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Abstract: Using soft-collinear effective theory (SCET), which provides a unified frame- 
i-C work for factorization, resummation of logarithms, and incorporation of universal nonper- 

turbative functions in hard-scattering QCD cross-sections, we present a new prediction of 
angularity distributions in e + e~ annihilation. Angularities r a are an infinite class of event 
& shapes which vary in their sensitivity to the substructure of jets in the final state, controlled 

by a continuous parameter a < 2. We calculate angularity distributions for all a < 1 to 
first order in the strong coupling a s and resum large logarithms in these distributions to 
next-to- leading logarithmic (NLL) accuracy. Our expressions for the next-to-leading order 
(NLO) 0(a s ) partonic jet and soft functions in the factorization theorem for angularity 
ff") distributions are given for the first time. We employ a model for the nonperturbative soft 

t-h function with a gap parameter which cancels the renormalon ambiguity in the partonic 

soft function. We explore the relation between the SCET approach to resummation and 
<^> past approaches in QCD, and discuss the advantages of the effective theory approach. In 

addition, we draw from the NLO calculations of the jet and soft functions an intuitive 
lesson about how factorization breaks down in the effective theory as a — ► 1. 
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1. Introduction 



Event shapes probe the hadronic final states produced in hard scattering processes for jet- 
like structure [1]. Two-jet event shapes e in hadronic e + e~ annihilations are constructed so 
that one of the kinematic endpoints corresponds to the limit of two back-to-back perfectly 
collimated jets. Different event shapes vary in their sensitivity to particles close to or far 
away from the jet axis and thus used in tandem probe the substructure of jets [2, 3]. Some 
examples of two-jet event shapes are the familiar thrust [4, 5], jet masses [6, 7, 8], and 
jet broadening [9], and the more recently introduced angularities [10]. The shape of the 
distributions in these variables depend on several energy scales, namely, the scale Q of the 
hard scattering, the scale of the invariant mass or typical transverse momentum of the 
jet hj, and the scale Aqcd of soft radiation from the jets involved in color recombination 
occurring during hadronization. Event shapes thus probe the behavior of QCD over a 
large range of energy scales, and indeed have been the source of some of the most precise 
extractions of the strong coupling constant a s [11, 12, 13]. 

Dependence on strong interactions at soft scales near Aqcd where QCD is nonpertur- 
bative would render predictive calculations impossible, without the use of factorization. 
Factorization separates an observable into pieces depending on each individual relevant en- 
ergy scale. Those pieces depending on large scales can be calculated perturbatively, while 
those depending on soft scales remain nonperturbative. If these soft functions are, however, 
universal among different observables or physical processes, then calculations of the factor- 
ized observables become predictive. A large number of two-jet event shape distributions in 
e + e~ annihilation can be factorized into hard, jet, and soft functions: 

— = H(Q;^) / dei de 2 de s Ji(ei;/i) J 2 (e25 Ai)5*(e s ;^)(5(e - ei - e 2 - e s ) , (1.1) 

otot de J 

where H(Q; (i) is the hard coefficient dependent only on the hard scattering e + e~ — > qq 
at center-of-mass energy Q, J\p are jet functions describing the perturbative evolution 
of the initially produced partons q, q into collimated jets of lower-energy partons, and 
finally S(e s ; fi) is the soft function describing the color exchange between the two jets 
leading to the hadronization of their constituent partons. This description introduces 
dependence on a factorization scale fi, at which the cross-section is factorized, into each 
of the individual functions. This dependence must cancel in the whole combination in 
Eq. (1.1). The full distribution and the individual jet and soft functions contain terms of 
the form (l/e)a™ ln m e which become large in the two-jet limit e — ► 0. The dependence 
of the hard, jet, and soft functions on the factorization scale \x can be determined from 
renormalization group equations, which can be used to resum the large logarithms [14]. 

The formidable achievements of proofs of factorization theorems for hard scattering 
cross-sections in QCD span a long and monumental history [15, 16]. More recently many of 
these theorems were reformulated in the language of soft-collinear effective theory (SCET) 
[17, 18, 19, 20]. This was done for two-jet event shapes for light quark jets in the series 
of papers [21, 22, 23, 24] and for top quark jets in [25, 26]. Some of the relations between 
the full and effective theory formulations of factorization were explored in [23, 24]. Equiv- 
alent results can be formulated in either language, although our discussion below will be in 
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the context of SCET, which we find advantageous for its intuitive framework for separat- 
ing physics at hard, collinear, and soft scales and its explicit Lagrangian for interactions 
between collinear and soft modes. These features facilitate the implementation of factor- 
ization and resummation of logarithms of ratios of all the relevant energy scales. At the 
same time that the effective theory provides us an intuitive framework in which to analyze 
the behavior of event shape distributions, the properties of the angularities themselves will 
in turn illuminate properties of the effective theory, and in particular, the conditions under 
which it is valid for the observables under consideration. 

To describe the conditions under which the distribution in a particular event shape 
factorizes as in Eq. (1.1), it is useful to write event shapes in a generic form. Many event 
shapes can be written in the form, 

e(X) = ^\p i L \f e ( Vi ), (1.2) 

where the sum is over all particles i in the final state X, p l ± is the transverse momentum 
of the ith particle and rji its rapidity relative to the thrust axis. Each choice of the weight 
function f e determines a different event shape. For example, for the thrust and jet broad- 
ening, /i_t(t/) = e"'^ and /b(t/) = 1. A continuous set of event shapes which generalize 
the thrust and jet broadening are the angularities r a [10], corresponding to the choice 

UM = z- M(l - a \ (1-3) 

where a is any real number a < 2. For a > 2, the function in Eq. (1.3) weights particles 
collinear to the thrust axis too strongly and makes the quantity Eq. (1.2) sensitive to 
collinear splitting, and thus not infrared-safe. The factorization theorem Eq. (1.1), however, 
is valid only for a < 1. At a = 1, the distribution of events in t\ is dominated by 
jets with invariant mass of order Aqcd- Thus, the jet and soft scales coincide, and the 
distribution cannot be divided into separately infrared-safe jet and soft functions, at least 
in the traditional form of the factorization theorem. This breakdown can be seen in the 
uncontrollable growth of a number of nonperturbative power corrections as a — > 1 [10, 23], 
or in the failure to cancel infrared divergences in the perturbative calculation of the jet 
or soft functions in the same limit, as we have recently explored in Ref. [27]. We review 
this breakdown of factorization in the explicit perturbative calculations we perform below. 
Any choice of weight function f e that sets a jet scale at or lower than the soft scale will 
ruin the factorization Eq. (1.1). 

The distributions for which the factorization in Eq. (1.1) breaks down might still 
factorize in a different form, by distinguishing collinear and soft modes not by their invariant 
mass, but by their rapidity, as proposed in [28]. We do not, however, pursue such a strategy 
here, and focus only on angularities with strictly a < l. 1 

The soft function evaluated at a scale fx s ~ Aqcd is nonperturbative. Evaluated at a 
higher scale, however, it can be calculated in perturbation theory. An appropriate model 

1 Even though traditional factorization breaks down for a — 1 (jet broadening), the resummation of jet 
broadening in QCD was performed in [9, 29] and nonperturbative effects were discussed in [30]. 
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for the soft function should interpolate between these two regimes. In our analysis we adopt 
a model like that proposed for hemisphere jet masses in [31] and for 6-quark distributions 
in [32], in which the soft function is a convolution, 

S(e s ; M ) = J de' s S PT (e s - e' s ; M )/ exp (e's " A e ) , (1.4) 

where 5" PT is the partonic soft function calculated in perturbation theory, and / exp is a 
nonperturbative model function. The gap parameter A e , proposed in Ref. [31], enters / exp 
through a theta function 9{e s — A e ) so that the minimum possible value of an event shape 
e of final states is A e , which is zero in the partonic distribution, but is nonzero due to 
hadronization in the actual distribution. The full soft function S(e s ;n) inherits its scale 
dependence from 5 PT (e s ; fi) and thus has a well-defined running with the scale \i. 

The partonic soft function S' PT (e s ; fi) contains a renormalon ambiguity due to the be- 
havior of the perturbative series at high orders. This ambiguity should not be present in 
the full physical distribution or the soft function, so the ambiguity in S PT is canceled by 
a corresponding ambiguity in A e . Shifting from A e to a renormalon- free gap parameter 
A e (/i) = A e — <5 e (/i) removes the ambiguity from the entire soft function Eq. (1.4). This 
greatly reduces the uncertainty in the predicted distribution due to such renormalon ambi- 
guities. These features were demonstrated in [31] for jet mass and thrust distributions. In 
this paper, we extend the soft function model and demonstrate that a similar cancellation 
occurs for angularities r a . 

Many studies of nonperturbative soft power corrections in event shape distributions 
have been based on the behavior of the perturbative expansions of the distributions, either 
the behavior of their renormalon ambiguities [33, 34] or their dependence on a postulated 
"infrared" effective coupling a s at low scales [35, 36, 37]. In particular, they led to the 
proposal of a universal soft power correction to the mean values of event shape distributions 
in the form [36, 37] 

(e) = <e>PT + "g 4 , (1-5) 

where (e)px is the mean value of the partonic distribution, and the coefficient of the 1/Q 
power correction is an exactly-calculable number c e dependent on the choice of event shape 
multiplied by an unknown nonperturbative parameter A which is universal for numerous 
event shape distributions. In [23] the operator definition of the soft function in the fac- 
torization theorem Eq. (1.1) was used to prove the relation Eq. (1.5) to all orders in a s . 
For angularities, c Ta = 2/(1 — a). This scaling of the power correction with a was observed 
in [38] based on the behavior of the resummed perturbative series for angularity distribu- 
tions after imposing an IR cutoff on the scale in a s (fi) and in [39] based on analysis of the 
distributions using dressed gluon exponentiation [40]. Below we will review the proof of 
the scaling in [23] based on the operator definition of the soft function independently of 
its perturbative expansion, and later use the scaling rule to constrain the nonperturbative 
model we adopt for the soft function in angularity distributions. 

The history of calculating event shape distributions using perturbation theory in QCD 
goes all the way back to QCD's earliest years. The thrust distribution for light quark jets 



- 4 - 



to 0(a s ) was calculated in [41], to which our fixed-order results for da/dr a reduce at a = 0. 
The resummation of the thrust distribution to NLL was performed in QCD in [42, 43] and 
to LL in SCET in [44, 45] (and later extended to N 3 LL in [11]). Our results are consistent 
with these SCET results at the appropriate orders for a = 0. The jet mass distribution 
for top quark jets was calculated and resummed to the same order in [26], with which we 
agree on the SCET jet and soft functions for a = in the limit mt = 0. The jet and soft 
functions for thrust or jet mass distributions can be derived easily from the "ordinary" 
SCET jet function J(k + ), and the hemisphere soft function S(k + , k~), because the thrust 
and jet mass depend only on a single light-cone component of the total four-momentum in 
each hemisphere (cf. [42]). These standard jet and soft functions were calculated to two- 
loop order in [46, 47, 48]. Angularities for arbitrary o, however, depend on both light-cone 
components k^ in each hemisphere, thus requiring the new calculations we perform below. 

In the original introduction of the angularities r a [10] the resummation of logarithms 
was achieved to the same next-to-leading-logarithmic (NLL) order that we achieve below, 
but without full inclusion of next-to-leading-order (NLO) jet and soft functions for the 
r a -distribution, which we calculate explicitly here for the first time. This improves the 
accuracy of our result for small r a . Our result is also improved in this region by adopting 
the soft function model Eq. (1.4) which cures unphysical behavior of the point-by-point 
distribution d<r/dr a as r a — > due to renormalon ambiguities. The results of [10] converted 
to the traditional form of an NLL resummed event shape distribution [42] were subsequently 
matched to fixed-order QCD at 0{a 2 s ) numerically in [38], improving the accuracy of the 
large- r a region. We perform this fixed-order matching only at 0(al). 

Comparing our result to those of [10, 38] elucidates the relation between SCET and 
traditional QCD-based approaches to resumming logarithms more generally. While the ad- 
vantages of SCET in achieving factorization or resummation of logarithms through renor- 
malization group evolution can of course be formulated without the explicit language of the 
effective theory (see, e.g., [10, 14]), the effective theory nevertheless unifies these concepts 
and methods in an intuitive framework that, we have found, allows us greater facility in 
improving the precision and reliability of our predictions of event shape distributions. Even 
though we do not go beyond the existing NLL resummation of logarithms of r a [10, 38] , the 
flexibility in the effective theory to vary the scales /j,h,j,s, where logarithms in the hard, 
jet, and soft functions are small and from which we run each function to the factorization 
scale fj,, allows additional improvements. For example, we are able to avoid any spurious 
Landau pole singularities which the traditional approaches usually encounter. (For previ- 
ous discussions on how the effective theory avoids spurious Landau poles present in the 
traditional approach, see Refs. [49, 50, 51].) 

The plan of the paper is as follows. In Sec. 2, we review the demonstration of factor- 
ization of event shape distributions in the formalism of SCET that was presented in [24] , 
recalling the introduction of the event shape operator e that returns the value of an event 
shape e of a final state X, constructed from the energy-momentum tensor. In Sec. 3, we 
calculate the jet and soft functions appearing in the factorization theorem for angularity 
distributions for a < 1 to one-loop order in a s . We recall the observations of [27] about 
how the breakdown of factorization as a — > 1 is observed in the infrared behavior of these 
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functions in perturbation theory. In Sec. 4 we solve the renormalization group equations 
obeyed by the hard, jet, and soft functions and resum leading and next-to-leading loga- 
rithms of r a in the perturbative expansions of these functions, and explain how we match 
the resummed distributions onto the fixed-order prediction of QCD at 0{a s ). In Sec. 5 we 
construct a model for the soft function in angularity distributions for all a < 1, based on 
existing models for hemisphere and thrust soft functions which contain a nonperturbative 
gap parameter introduced in [31], which cancels the renormalon ambiguity in the pertur- 
bative series for the soft function. In Sec. 6 we present plots of our final predictions of 
angularity distributions using all the results of Sees. 3-5. In Sec. 7 we compare and con- 
trast the SCET approach to predicting resummed angularity distributions to those based 
on factorization and RG evolution in full QCD [10] and to the traditional approach to re- 
summation [38, 42]. In Sec. 8 we present our conclusions, and in the Appendices, we verify 
a consistency relation among the hard, jet, and soft anomalous dimensions for arbitrary 
a, provide some technical details necessary for the solution of the RG equations for the 
jet and soft functions, and explain our procedure to calculate angularity distributions at 
fixed-order in QCD at 0(a s ), noting the hitherto unnoticed property of the angularities 
that they fail to separate two- and three-jet-like events for values of a < —2, and so cease 
to behave exactly as "two-jet" event shapes. 

2. Review of Factorization of Event Shape Distributions 

We begin by reviewing the factorization of event shape distributions in the formalism of 
SCET, presented in [24]. 

2.1 Event shape distributions in full QCD 

The full QCD distribution of events in e + e~ — ► hadrons in an event shape variable e is 
given, to leading-order in electroweak couplings, by 



H7 = ^E M 4 ^' E Li lu (0\jf(x)\X)(X\jr(0)\0)S(e-e(X)), (2.1) 



where q = (Q, 0) is the total four-momentum in the center-of-mass frame, the sum is over 
final states X, and e{X) is the value of the event shape e of the state X. The final state 
is produced by the vector and axial currents, 



where Ty = 7^ and = 7^7 5 and the sum is over quark flavors / and colors a. The 
leptonic tensor, which includes contributions from an intermediate photon and Z boson, is 
given by 




i=V,A 




(2.2) 




Q 2 



) 




2Q 2 v e v f Q f Q> e 2 + a 2 e )v 2 f 
QZ-Ml + {Q 2 -MlY 



(2.3a) 




(2.3b) 
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where Qf is the electric charge of / in units of e, and Vf,a,f are the vector and axial charges 
of/, 

v f = ^-fT ^(7/ -2Q /S in 2 M, a f = ^-lT 7T T r ^ 

2 sm By/ cos Ow 2 sin Vyy cos 6w 

As shown in [24], the sum over hadronic final states remaining in Eq. (2.1) can be performed 

by introducing an operator e that gives the event shape e{X) of a final state X. This 

operator can be constructed from a momentum flow operator, which in turn is constructed 

from the energy-momentum tensor. That is, 

e \X) = e(X) \X) = i f°° dr, f e (v)£ T (v, *) I*) , (2-5) 

J — oo 

where t is the operator yielding the thrust axis of final state X, and £t(V'i Q is the transverse 
momentum flow operator, yielding the total transverse momentum flow in the direction 
given by rapidity 77, measured with respect to the thrust axis, in a final state X, 

£ T (v,i) \X) = / d<p lim R 2 / dthiT 0l (t,Rh) \X) = V \p l ± \5( V - m ) \X) , 

cosh° i] J R^oc J f^x 

(2.6) 

which is closely related to the energy flow operator proposed in [52]. The thrust axis 
operator t can be constructed explicitly, as shown in [24]. After matching onto SCET, 
however, an explicit construction is not necessary, as the thrust axis is simply given by 
the jet axis n labeling the two-jet current. The difference between the two axes introduces 
power corrections in A which are subleading, as long as a < 1 [10, 24]. Using the operator 
e, we perform the sum over X in Eq. (2.1), leaving 

!r = ^2 / d ' xeiq ' X E L ^ <°l ^VWe - e)jT(O) |0) • (2.7) 
ae 1(4 J i=V,A 

2.2 Factorization of event shape distributions in SCET 

To proceed to a factorized form of the distribution Eq. (2.7), we match the current and 
the operator e onto operators in SCET. To reproduce the endpoint region of the two-jet 
event shape distribution, we match the QCD currents $ onto SCET operators containing 
fields in just two back-to-back collinear directions, 

3i( x ) = E E Cnn{Pn,Pn] Ll)O n n{x]p n ,Pn) , (2.8) 
n Pn,Pn 

summing over the direction n of the light-cone vectors n,n = (1, ±n), and label momenta 
PrnPn- The two-jet operators [22, 53], after the BPS field redefinition [20] with soft Wilson 
lines, are 

O nii (x;p n ,Pn) = ^'^Xn,^ {x)Y n (x)T? Y n ( X ) X n, Pfi (x) , (2.9) 

where Ty = 7^ and = 7^75. The soft Wilson lines are the path-ordered exponentials 
of soft gluons, 



Y n (x) = Pexp 



/*oo r /'oo 

ig \ n ■ A s {ns + x) , Yn(x) = Pexp ig / n ■ A s (ns + x 
Jo J I Jo 



(2.10) 
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with A S ,A S respectively in the fundamental or anti- fundamental representation. The jet 
fields Xn = Wtt, n and Xn = W ] n i n are combinations of collinear quark fields made invariant 
under collinear gauge transformations by Wilson lines of collinear gluons [18, 19], where 



W n (x) = ^2 ex P 



perms 



(2.11) 



where q is the label momentum of the collinear gluon field A n , and V is a label momentum 
operator which acts as VA n ^ q = (n • q)A nA [19]. Recall that, in SCET, collinear momenta 
Pc = P^ + k^ are divided into a large label piece, p^ = (n • p)n^ /2 + p^_, and a residual 
piece, k 11 , where n ■ p is 0(Q), p± is 0(QX), and k is 0(QA 2 ). The residual momenta are 
the same size as soft momenta, k s , of 0(QX 2 ). Below, however, we will see how the natural 
scaling of the collinear modes varies with the choice of observable r a . The integral over x 
in Eq. (2.7) enforces that the label momenta of the jet fields in the two-jet operator satisfy 
n-p n = -n-pn = Q and p^ = p± = 0. 

We must also match the operator e in full QCD onto SCET. To do so we simply 
replace the QCD energy- momentum tensor appearing in the definition Eq. (2.6) with 
the energy-momentum tensor in SCET, and, as noted above, set the thrust axis equal to 
the jet axis n in the two-jet operator O n n- After the BPS field redefinition, to leading order 
in A the SCET energy-momentum tensor is a direct sum over contributions from fields in 
the n, n collinear and soft sectors, since the Lagrangian splits into these separate sectors 
with no interactions between them. (Beyond leading order in A, there are power-suppressed 
terms in the SCET Lagrangian in which interactions between collinear and soft fields do 
not decouple following the BPS field redefintion [54, 55, 56, 57].) Then the event shape 
operator e splits into separate collinear and soft operators, 

e = e n + e n + e s , (2.12) 

where each ei is constructed only from the energy-momentum tensor of sector i of the 
effective theory. So, finally, the event shape distribution in SCET factorizes into purely 
hard, collinear and soft functions, 

1 do - f 

— — = H(Q;n) / de n de n de s 5(e - e n - e n - e s )J n (e n ;iJ,) Jn{e n ;^)S(e s ;fi) , (2.13) 
(To de J 

where the hard coefficient is the squared amplitude of the two-jet matching coefficient, 

H(Q-n) = \C nn (Qn/2,-Qn/2;ti)\ 2 , (2.14) 

and the jet and soft functions are given by the matrix elements of collinear and soft oper- 
ators, 

S(e s ; M ) = ^ Tr <0| Y{(0)Y^0)5(e s - e s )Y n (0)Y n (0) |0) , (2.15) 

and 

/dl + f dk~ 

-^—Jn(e n ,l + ;fJ>), Jn(en,fJ>)= / -7; — Jni&n, k~ J £t) , ( 2 -16) 
27T I ITT 
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where 

J n (e n , = ~rf- Tr / d 4 x e ilx (0| Xn, Q (x) a S(e n - e n ) X n, Q (0)p |0) (2.17a) 



Jfl(eB,A;-;M) (f) = jVc ^ / ^^' X (0\Xn,-Q(x)p5(e n - e n )Xn,- Q (0) a \0) . (2.17b) 

In Eqs. (2.15), (2.17a), and (2.17b), the traces are over colors. Also, in Eq. (2.13), we have 
divided the distribution by the total Born cross-section for e + e~ — ► qq, 



4na 2 N c 

£J = 



E 



3Q , 



2 2Q 2 v e v f Q f Q\vj + ai)(vj + aj) 
f Q*-M* + (Q 2 -Af|) 2 



(2.18) 



The n-collinear jet function J n depends only on the l + = n ■ I component of the residual 
momentum, and J n on k~ = n ■ k, as only the n ■ d derivative appears in the n-collinear 
Lagrangian, and n-d in the n-collinear Lagrangian, at leading order in A [18]. In angularity 
distributions, the jet functions are independent of the residual transverse momenta k±, l± 
as long as a < 1 [24]. 

In Sees. 3 and 4 we calculate the above hard, jet, and soft functions for angularity 
distributions to next-to-leading order in a s , and solve for their dependence on \x through 
the renormalization group equations, which will allow us to sum large logarithms of r a . 

2.3 Universal first moment of the soft function 

As shown in [23], the behavior of the soft function Eq. (2.15) under Lorentz boosts in the 
n direction implies a universal form for its first moment. The vacuum |0) and the Wilson 
lines Y n ^n(0),Yn t n(0) are all invariant under such boosts, while the transverse momentum 
flow operator £t(v) appearing in the definition of e s transforms as £t(v) ~* £t(v') under a 
boost by rapidity rj — rj. These properties imply that the first moment of S(e s ; n) is given 
by 

de s e s S(e s ;u) = ^M, (2.19) 



/ 



dvfeiv) (2-20) 



where 

1 

c P = — 

Q 

Afi) = ^ Tr <0| H(0)^(0)£r(0)rn(0)F fl (0) |0) . (2.21) 

The coefficient c e is exactly calculable from the definition of the event shape e in Eq. (1.2) 
while A(n) is not fully calculable due to the contribution of nonperturbative effects, but is 
completely independent of the choice of variable e. The first moment Eq. (2.19) is universal 
for all event shapes of the form Eq. (1.2) in this sense. For angularities, using Eq. (1.3) 
and Eq. (2.20), 

poo o 

drje -\r)\(l-a) = ( 222 ) 

oo 1 - a 

This scaling of the first moment of the soft function for angularities will constrain the 

parameterization of the nonperturbative model for the soft function that we introduce in 

Sec. 5. 
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3. Fixed-order Perturbative Calculations of Hard, Jet, and Soft Functions 



In this section we calculate at next-to-leading order, that is, 0(a s ), in perturbation theory 
the hard, jet, and soft functions, H(Q; fj,), J^ n {ja' n \ /u), and S a {r^\ ^), in the factorization 
theorem for angularity distributions, which is given by Eq. (2.13) with e = r a . 2 

3.1 Hard function at NLO 

The hard function H(Q;fj,), given by Eq. (2.14), is the squared amplitude of the two-jet 
matching coefficient C n n(Q, —Q', [i)- This matching coefficient was calculated, for example, 
in [49] in the context of DIS and in [22] for e + e~ annihilation, to NLO. It is found by 
calculating a matrix element of the QCD current Eq. (2.2) and SCET current Eq. (2.8) 
(for example, {q(p q )q(pq)\ |0)), and requiring that the two match. Since the matching of 
the currents is independent of the observable being calculated, we do not need to repeat 
the matching calculation here, and simply quote the result. The matching coefficient 
Cnn(Pn,Pn, fi) in the SCET current Eq. (2.8) is given by 



CnniPmPn] /^) — 1 ~. 

47T 



7T 2 , 2 / fj 2 \ , / I? 

+ ln 2 _ - +31n ' P 



6 \2pn-pnJ \2p n -p n 



(3.1) 



Here and in the remainder of this section, a s = a s (/z). The hard function H(Q;fi) in 
Eq. (2.14) is thus 

cx s Cf ( 7tt 2 , 9 u 2 , u 2 



H(Q;,,) = l--^(8-— + ln 2 ^+31n^) . (3.2) 



The additional contribution to the coefficient of 7r 2 in going from Eq. (3.1) to Eq. (3.2) is 
due to the sign of 2p n ■ p n = —Q 2 , following the conventions of [19]. 

The bare SCET two-jet operators in Eq. (2.9) are renormalized by the relation 

Onn( x 'iPn,Pn) = Z (p n , Pn\ V) O nfl {x] p n , pn) , (3.3) 

where the renormalization constant, calculated using dimensional regularization to regulate 
the UV divergences in d = 4 — 2e dimensions, is given by 



Zoipn,Pn\ H) = 1 + 



2 , 2, / \ 3 



e 2 + e ln \2p n - pn) + e 



(3.4) 



Matching the QCD current Eq. (2.2) onto only two-jet operators in SCET is sufficient to 
describe accurately the two-jet region near r a = of angularity distributions. To calculate 
accurately also the tail region to 0(a s ), where the jets broaden and an additional jet begins 
to form, we would need to include a basis of three-jet operators in Eq. (2.8) as well [44, 58]. 
But since we are mainly interested in obtaining the correct shape of the two-jet region, 
we do not pursue this approach here. We will simply calculate the whole distribution in 
SCET with only two-jet operators, and then match the tail region numerically onto the 
fixed-order prediction of full QCD. This will be described more precisely in Sec. 4.4. 



2 Note that here and below a superscript n on a quantity is not a power but denotes "n-collinear" just 
as n denotes "n-collinear" and s denotes "soft". 
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3.2 Cutting rules for weighted matrix elements 

The jet and soft functions that typically appear in factorizations of hard cross-sections in 
SCET are defined in terms of matrix elements of the products of collinear and soft fields, 
which are related to the imaginary part of the matrix element of a time-ordered product 
of the fields according to the optical theorem, 

d 4 se i «- <B <O|0(a;)0 t (O)|O) = DiBc J d 4 x e^O |T ^(x)^ f (0)|0> . (3.5) 

The right-hand side is then related to the sum of all cuts of the relevant Feynman diagrams 
using the standard Cutkosky cutting rules. 

However, for more generic jet observables such as angularities for the jet and soft 

functions that appear in factorization proofs contain matrix elements in which additional 
operators are inserted between the collinear and soft fields in the definition of the traditional 
jet and soft functions [59]. For the matrix elements involving the extra insertion of such 
operators, we need to generalize the cutting rules for calculating these matrix elements 
from Feynman diagrams. 

For the case of angularities, the jet and soft functions given in Eqs. (2.15), (2.17a), and 
(2.17b) differ from the traditional jet and soft functions by the insertion of the delta function 
operator 5(r a — f a ). We denote the appropriate generalized prescription for calculating the 
new matrix element from the Feynman diagrams of time-ordered perturbation theory as 
the "r a -discontinuity," 

d 4 xe^ x (0\<f>(x)S(T a -f o )<£t(0)|0) = Disc TQ /Vxe^(O|T<^(x)0t(O)|O) . (3.6) 

The Disc Ta prescription is to cut the diagrams contributing to the matrix element of time- 
ordered operators just as for the usual matrix elements in Eq. (3.5) but to insert an ad- 
ditional factor of 5{r a — T a (X)) for each cut, where X is the final state created by the 
cut. 3 This prescription corresponds to reinserting a sum over a complete set of final states 
between the delta function operator and <^(0) in Eq. (3.6), and is precisely how we would 
calculate the full differential cross-section as written in Eq. (2.1). In the next two subsec- 
tions we illustrate extensively the use of the Disc Ta prescription. 

3.3 Calculation of the soft function to NLO 

The diagrams that contribute to the soft function are shown in Fig. 1. From Eqs. (1.3), 
(2.5), and (2.6), the contribution to the angularity from an on-shell soft gluon with mo- 
mentum k is 

T .jM e -^\j f - "->- k \ ( 3.7) 

Q I ilArl ^Ifc+I* f or u + >k- 



3 The operator-based method that was developed in [60] for calculating weighted cross-sections can be 
used to relate matrix elements such as in the left-hand side of Eq. (3.6) directly to the ordinary discontinuity 
of matrix elements of time-ordered products of fields. However, for the scope of this paper, we choose simply 
to apply the prescription Eq. (3.6). 
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(A) (B) (C) (D) 

Figure 1: The (A), (B) real and (C), (D) virtual contributions to the soft function. The gluons 
all have momentum k. 



Since cutting a gluon puts it on shell, the operator returns these values when acting on 
a cut soft gluon. When no gluon is in the final state cut, the operator simply returns 
zero. The real and virtual diagrams then contain delta functions, which we denote 5r and 
5y, respectively, 

6r = S r (t°, k) = 0(AT - k + ) 8(r s a - i \k+\ 1 ^ \k-\%) (3.8a) 

+ e(k + -k-)5(r°-^\k-\ 1 -^\k + \^, 
S V = 5 v {r s a ) = 5{r s a ) . (3.8b) 
In terms of these delta functions, the (bare) perturbative soft function can be written 




Sr + 2 




(3.9) 



where we used that the tree-level contribution is just 8(t%) and that the two real and the 
two virtual diagrams in Fig. 1 give identical contributions. 

In pure dimensional regularization, the virtual contributions are scaleless and hence 
vanish so we only need to evaluate the real diagrams. They add to 




>S R = 2g 2 f j 2e C F n-n 
Performing the k integrals gives 
S p a T{0 \T s a ,n)=5(T°) + 9(T l 



d d k 1 1 
(2ir) d A^AF 



2^5{k~k + - \k ± \ 2 )6(k-)5 R (T°,k) . (3.10) 



a s Cpn-n f Ait^ 



2\ e 



1 1/1 



l+2e 



(3.11) 



vr(l-a) \ Q 2 J r(l-e)e 

Nonzero values of regulate the IR divergences, and so here the 1/e pole is of UV origin, 
e = euv- 

Applying the distribution relation (valid for e < 0) 



0{x) 5{x) 



x 



l+2e 



2e 



+ 



'00*0 



2e 



(x) In x 



x 



+ 0(e 2 ) , 



(3.12) 
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where 



'e(x)\n n ( x y 


= lim 


X 





x n + 1 



to Eq. (3.11) we obtain the final result for the (bare) angularity soft function, 



(3.13) 



(3.14) 



where to NLO the renormalized soft function, S% , is given by 



1 



a s C F ( 1 ^ ^ 2 



7T 



tt(1 - a) \2 Q 2 12 
and the renormalization factor, Zs, is given by 



+ 



2a s C F 
7r(l - a) 



Z s (r-^) = 5(r a s ) 



_ a s C F / 1_ 1 /x 2 
vr(l-a)Ve 2 e Q 2 



+ 



1 2a s C F 
e 7r(l — a) 



(3.15) 



(3.16) 



3.4 IR structure of the soft function 



While the mathematical identity in Eq. (3.12) allowed us to arrive at our final result, 
Eq. (3.15), the origin of the 1/e poles became obscured through its use. In fact, the use of 
Eq. (3.12) is only valid for e < which suggests that the 1/e pole on the right-hand side of 
Eq. (3.12) is of IR origin. The virtual diagrams, while formally zero in pure dimensional 
regularization, play the role of converting this IR divergence into a UV divergence by adding 
a quantity proportional to (1/euv — V e iR.) to the coefficient of 5(r*), if the final result is in 
fact free of IR divergences. Naively it seems that this conversion cannot possibly occur for 
arbitrary a, because the 1/e poles in the real diagrams have a-dependent coefficients (see 
Eq. (3.16)), while the virtual diagrams contain no apparent a dependence. Nevertheless, 
by carefully examining the contribution of both the real and virtual diagrams, we will 
show that, for a < 1, the virtual diagrams play precisely this role and convert each IR 
divergence in the real graphs into UV, but that for a > 1, this cancellation is incomplete. 
This is accomplished through an analysis of integration regions in the loop momentum 
integrals that avoids the use of explicit IR regulators. Our presentation here complements 
our discussion of these issues in [27]. 

Using that J^dx [ln n (x)/x] + = 0, the contribution to the coefficient of 6(r^) can be 
isolated by integrating the diagrams over from to 1. We find that the contribution 
from the real diagrams can be written as 



dr* 




a s CFn-n (47r/x 2 )' 



dk + dk~{k + k~) 



-l-e 



2vr T(l-ejjn 

where 1Z is given by the region of positive k + and k~ such that 

(AT) § (A;+) 1 -f < Q for AT > k + 
(A; + )f (A;-) 1 ^! < Q for AT < k + . 



(3.17) 



(3.18) 
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Figure 2: The regions of integration for the coefficent of <5(r^) in Sa ( T a) m the (A), (B), (C) 
k~,k + and (D) fc~, planes. The regions of integration for both (A) the real contribution TZ and 
(B) the virtual contribution V contain both UV and IR divergences. Since the integrands for the 
two contributions differ only by an overall minus sign, (C) the region resulting in their sum S, is the 
complement of TZ and contains only UV divergences for a < 1. The dashed line in (C) represents 
the line of constant k + k~ = Q 2 . 



This region is plotted in Fig. 2A for various values of a. 

The contribution of the virtual diagrams to the coefficient of 5(t^) sums to 



dr„ s 




, c\ d k 

2g 2 ^C F n-n ' 



(2vr) d k~ - i0+ k+ + i0+ k+k~ - |kj 



2tt 



T(l - e) 



dk + dk~(k + k 



-\-l-e 



- i0 + 

(3.19) 



where V is the entire positive k + , k~ quadrant, plotted in Fig. 2B. 

The two contributions to 8(t^), Eqs. (3.17) and (3.19), are each both UV and IR 
divergent, but as we will show, their sum is convergent for e > and so is only UV 
divergent. Since the form of the integrand is the same and the virtual contribution differs 
only by an overall minus sign, it converts the region of integration of the real contribution, 
TZ, into the complementary part of the positive k + ,k~ quadrant (see Fig. 2) which does 
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not include the IR divergent regions fc+ -> 0. Note that as a — ► 1, the boundary of the 
region of integration 1Z approaches the curve of constant k + k~ = Q 2 . With this boundary, 
the integral over the region S does not converge for either positive or negative e, implying 
that both IR and UV divergences are present. 

That the region S has only UV divergence for a < 1 and has both UV and IR divergence 
for a = 1 is perhaps more clearly seen in the k~, plane. The integral of the soft diagrams 
over in terms of these variables is given by 



dr„ s 




Sr + 2 




a s CFri-n (47r^ 2 ) e 



2tt T(l-e) Js jfe- (k3)i+e ' 



dArdk 2 



and the resulting region S in terms of k and k^ for a < 1 is 



"2(l-«) 



< 



jfe- 

~Q 



< 



Q 2 



2-a 
2(l-a) 



with k 2 ± > Q 2 



The region S is plotted for several values of a in Fig. 2D. The limiting case a 
includes the IR divergent region k~ — ► for all k^ > Q 2 . 
Performing the integral over S we obtain 



(3.20) 



1 clearly 



dT a s S(°)(T^) = l 



a s Cpn-n ( 4tt/j, 
2vr(l - a) V~Q2 



1 



s 2 r(i 



(3.21] 



After expanding Eq. (3.21) in e, we find that the coefficient of 5(t*) in Eq. (3.15) is 
unchanged, except that for a < 1 all the 1/e poles are unambiguously of UV origin. 

A lesson from this analysis is that in pure dimensional regularization, the coefficient of 
(1/euv — IAir) i n a virtual diagram cannot be determined from the virtual diagram alone, 
but only together with the real diagram whose IR divergence it is supposed to cancel. The 
reason that the virtual subtraction can depend on a even though by itself it is independent 
of a is that the area of overlap between the integration regions of real and virtual diagrams 
depends on a. 



3.5 Calculation of the jet functions to NLO 

Now we proceed to calculate the jet functions given by Eqs. (2.16) and (2.17). The diagrams 
that contribute to J" are shown in Fig. 3, and the Feynman rules necessary to calculate 
these diagrams are found in [18]. The total momentum flowing through each diagram is 
Qn/2 + I, with the label component Qn/2 specified by the labels on the jet fields in the 
matrix elements in Eq. (2.17a), and I the residual momentum. The total momentum of the 
gluon in each loop is q, which has both label and residual components. All results for the 
anti-quark jet function J" can be found from those for the quark jet function J" with the 
replacement n <-> n and so we calculate explicitly only J". 

Cutting the diagrams in Fig. 3 in all possible places, we can cut through the gluon loops 
or through one of the individual quark propagators connected to a current. We naturally 



- 15 - 



-® ®— 



-o ®~ 




(A) 



(B) 



(C) 



(D) 



Figure 3: Diagrams contributing to the angularity jet function J7™(t™, l + ) with incoming momen- 
tum I — ttQ + and gluon momentum q: (A) Wilson line emission diagram and (B) its mirror; 
(C) sunset and (D) tadpole QCD-like diagrams. The contributions to the jet function J"(t") 
are given by the integrals of these diagrams over the + component of the incoming momentum, 

call these classes of cut diagrams "real" and "virtual" respectively. The real and virtual 
diagrams contain the delta functions, 

l r. 



S R = S R (T^q t l+ 
5 v = S v (r^l+) = Mr, 



+ (Q 



i a 



(3.22) 



which are obtained using Eq. (3.7). In this case we simply consider the contribution to r a 
from a final state with a single on-shell collinear quark of momentum I for Sy and from a 
final state consisting of an on-shell collinear gluon of momentum q together with an on-shell 
collinear quark of momentum I — q for 5r, and use that the ' — ' component of momentum 
is always larger than the '+' component for on-shell collinear particles. The momentum I 
flowing through the diagrams in Fig. 3 has a label component which is fixed to be Qn/2 
by the labels on the collinear fields in the matrix element in Eq. (2.17a). 

Before turning to evaluate the diagrams in Fig. 3, we first perform a few simplifications 
to facilitate the computation. First, we note that the Wilson line emission diagram, Fig. 3A, 
and its mirror, Fig. 3B, give identical contributions. Second, we employ the fact that the 
number and complexity of jet function diagrams needed in loop calculations is reduced by 
noticing that the QCD-like diagrams can be computed using ordinary QCD Feynman rules 
with appropriate insertions of the projection operators P n = ^/4 and Pn = f/itf/A [47, 61]. 
In particular, for our one-loop example we use that the sum of Fig. 3C and Fig. 3D reduces 
to 



+ 



P / 5STO \ P- 



(3.23) 



Next, we relate the r^-discontinuity to the ordinary discontinuity, 



Disc T ; 
= 2 



Disc 



* + 1 

+ Pn 



+ Pn, 



Pn 



2/""% | ; +2P n jA f n 



Pn 



S v + 



(3.24) 



Sr-Sv) , 



where in the third line the we used that the real diagrams induced by taking the disconti- 
nuity in the first term cancel the coefficient of 5y in the second term. 
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Now, since <5y(r™, Z + ) has no dependence on the loop momentum q, it factors out of the 
d d q integrand. This implies that, after adding the tree-level contribution to the one-loop 
r"-discontinuity in Eq. (3.24), we can write the NLO jet function as 



J^°\t2,1+;^ = 2it5(1+)6(t:)^ +Disc r . 



where Jn\l + ;/j>) is the standard jet function [20], 
J^/^-^Disc 



(3.25) 



Pn 



5r - 5y 



N c 



Jd 4 xe u - x Tr(0\Txn^(x) x„,q(0)|0) 

^\ - + Pn r iT\ Pn] + 0{al) , (3.26) 



= 2vr5(/ + )^ + Disc 

containing no additional operator insertions. Each term on the second line of Eq. (3.25) is 
then well-defined 4 and straightforwardly calculable. In fact, J n (l + ; //) has been calculated 
to two loops [47], and we expect that the techniques we employed above are the most 
practical way to extend our results to two loops. The additional term on the second line of 
Eq. (3.25) is a sum of real emission diagrams containing a difference of the delta functions 
<5r and 5y. Note that for the special case a = 0, 5v(t", l + ) = $r(t2, q, l + ) and this 
additional term vanishes, so J n = J2=o- This is why only the standard jet function is 
needed when a = 0. 

To find the angularity jet function J™ (r^; //), we must integrate Eq. (3.25) over l + as 
in Eq. (2.16), 



•C (0) (Cm) 



dZ+ 
2?r 



(3.27) 



By integrating the known one-loop expression for J n °\l + ;fi) (see, e.g., [62, 63]), we find 
that the contribution of the first term in Eq. (3.25) is 



^(»)( i+ ; rtiv = JW){l + ^ 



4 3 4 , u 2 , 2 u 
^ + - + -ln^ + 21n 2 



e e 

,2 



Q 2 



+ 31n^ + 7-7T 2 

Q 

- + 3 + 8 In — — ^., , 
e Q(r^)V(2-o) 



1 - a/2 

It is well known that all 1 /e poles in this expression are of UV origin. 



(3.28) 



4 By this we mean that had we evaluated the individual cut virtual QCD-like diagrams contained in 
the first line of Eq. (3.25) directly, we would have encountered the complication of cutting one lone quark 
propagator and thus putting the second lone, uncut quark propagator on shell also. 
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We find that the term involving the real QCD-like diagram in Eq. (3.25) is 



/ 



dl+ 
2ir L 



P, 



Pm 



5r - Sy 



d d q 
(2vr) Q 



(r 



x ((-2ttz)%V - |qj 2 )#(<r)) ^(d R (r^q,l+) - 6 v (r^l + ) 

x ((-2vri)5((Q - q-)(l+ - q+) - \q ± \ 2 )6(Q - q~) 
a s C F 7^ /4vr/i 2 \ e 2(1 - e) / 1 \ 1+ ^72 



2vr(2-o)2 V Q 2 J r(l-e) \t£ 



dx x 



{x a ~ l + (1 - x)^ 1 ) - (x(l - x)Y 



where we defined x = q /Q. This expression is finite as e — > 0. 
For the term involving the real Wilson line diagram, we find 



df 
2tt L 



g 2 y? e Cp n-n 



S R - 5 V 

dl+ 1 f d d q 1 



-2m)5{q + q- - |q±| 2 )%~ 



2 J 2vr 1+ J (2ir) d q- 
(Q - q-)((-2m)5((Q - q~)(l + - q + ) - \q±\ 2 )9(Q - q 



-Q^-27Ti)6{Q(l + -q + )) 
The piece with 5r can be written as 



(s R (r:, q ,l+)-S v (r:,l+) 



I 



dl+ 
2tt 



a s C F n-ni f^fJpY 1 



2vr(2 - a) 2 \ Q 2 J T(l - e) \r* 
1 dx 



o x 



(1-x) {x a ~ l + (1 - a;) " 1 ) i-« 



'dx 



(3.29) 



(3.30) 



(3.31) 



-x i-"/ 2 



o x 



and the piece with 5y is 



/ 



dj+ 
2vr 



Sv = 9«) 



1 \ 1+ T^72 

2vr(2 - a) 2 V Q 2 J T(l - e) 



a s CFn-nr/i /47r/i \ 1 



dx 



(1 — x) (x(l — x))~ 



'dx 



(3.32) 



The second term in brackets in each of Eqs. (3.30), (3.31), and (3.32) corresponds to the 
zero-bin subraction [28] needed to avoid the double counting of soft modes [23, 64, 65]. 
Note that from the expressions in both Eqs. (3.31) and (3.32), the zero-bin contributions 
are scaleless and hence formally zero. Their role is to convert the IR divergence {q~ — > 0) 
in each integrand into a UV divergence {q~ — > oo) for a < 1. After this subtraction, both 
of the integrals over x in brackets are convergent for e > 0. 



- 18 - 



Subtracting Eq. (3.32) from Eq. (3.31) and performing the integral over x we find that 



'dZ+ 
~2tT 



Disc 



+ Pn 



Pn 



a s CF 4 ( 47T/X 



2vr(2 -a)2\ Q 2 J T(l - e) e \t% 



5 R - S v 
1 \ 1+ T^7 



(3.33) 



2a 
1 - a 



+ e- 



2a(vr 2 - 9) 
3(2 - a) 



1 - a/2 



o 



'dx 1 x + x2 l 2 l n [(i _ + x i-»] + o( e 3j 



where the overall 1/e pole is of UV origin from the discussion above. 

Applying the relation Eq. (3.12) to Eq. (3.33) and adding the result to Eq. (3.28), we 
arrive at our final expression for the (bare) NLO angularity jet function, 



4 n(0) (C /x) = J dTf Zj{t: - t:>; /i) 



(3.34) 



where the renormalized jet function, J™, is 



J a n (C/i) = o(r a ") jl + 
o. s Cf 



CXsCf 

3 



1 - a/2 1 2 fi 2 3 . ^ , 



7T 



1 



+ 



In 



4 1 - a/2 I- a Q(r") 1 /(2-a) y ^ T( 



(3.35) 



where we defined 



1 {7- 13a/2 _ vr 2 3 - 5o + 9a 2 /4 
/(a) = I^/2 V 4^ 12 T^a 



'dx 1 X + x2 / 2 l n [(i _ x )i~* + s i-»] 



(3.36) 



and the Z-factor is given by 



Zj(t2;ij,) = 5(t ( 



1 + 



o. s Cf / 1 — a/2 / 1 1 



7T 



1 - a 



+ - In 



+ 



4e 



1 oi s Cf 
e 7r(l — a) 



(3.37) 



3.6 IR structure of the jet functions 

As we showed in Sec. 3.5, the 1/e pole in front of the plus-distribution corresponds to a UV 
divergence. However, as we discussed for the case of the soft function in Sec. 3.4, the use 
of Eq. (3.12) means that we can not immediately make the same claim for the poles in the 
coefficient of 5(r"). We now perform an analysis similar to that in Sec. 3.4 by integrating 
over < r™ < 1 to isolate this coefficient and study its divergent structure in the resulting 
q~, integration regions, complementing our discussion in [27]. 

The diagrams (C) and (D) in Fig. 3, being equivalent to diagrams in full QCD as noted 
above, are manifestly infrared-finite and do not need to be analyzed in further detail. The 
Wilson line graphs (A) and (B) potentially contain infrared divergences that we must 
identify more carefully. 
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If the jet function is infrared-safe, infrared divergences in virtual and real diagrams, 
with proper zero-bin subtractions taken, should cancel and leave purely UV divergent 
integrals. The contribution of the sum of the real and virtual Wilson line diagrams to the 
coefficient of <5(r") in the jet function fl (t™) is 



1 

2 I dr a 

'0 



d/+ 

~2~7T~ 



S v + 



(3.38) 



a s C F (4V) £ 
vr r(i - e) 



dq dqj 
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2 M- 



dq dqj 



Jo 



(qi) 1 - 



where the last integral is the zero-bin subtraction of the naive collinear integral in the first 
term. The naive integration region J is shown in Fig. 4 and is given by < q~ < Q and 



ql> <Q 



1 



(Q 



\l—a 



+ 



1 



\1— a 



z/2 



The zero-bin region J§ is given by q > and 

qi > [Q{€ 



\1— a] l-a/2 



(3.39) 



(3.40) 



The resulting integral for the total contribution of the zero-bin-subtracted Wilson line 
diagrams to the coefficient of <5(r") in the jet function is 





fdl + ' 


[ dT 'J 
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S v + 



a s C F (47r/i 2 ) £ 
vr r(l-e) 



dq dq^ 
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(qi) 1 - 
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where the region J resulting from combining J and J7q, with a relative minus sign in the 
integrands, is also shown in Fig. 4. 

The shape of the final integration region J in Fig. 4 demonstrates that the scaleless 
virtual and zero-bin integrals succeed in converting IR divergences in the real diagram 
contributions into UV divergences for all a < 1. The integral over J in Eq. (3.41) converges 
for e > if and only if a < 1. The result of performing this integration, after including 
the contributions of the QCD-like diagrams in Fig. 3C and D, agrees with the coefficient 
of 5(t™) that is obtained by (naively) using the relation Eq. (3.12) in Eq. (3.33). 

3.7 Infrared safety, factorizability, and the effective theory 

In the one-loop calculations of soft and jet functions above, we observed that infrared 
safety of these functions, and, thus, factorizability of the angularity distributions, required 
a < 1. By analyzing explicitly the regions of integration over loop momenta in real and 
virtual graphs, we were able to identify when the loop integrals contained infrared or 
ultraviolet divergences. Cancellations of regions in real gluon diagrams sensitive to IR 
divergences relied crucially not only on the addition of virtual diagrams but also on zero- 
bin subtractions from collinear diagrams (see also examples in [28, 64, 65, 66]). 
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Figure 4: Regions of integration for the coefficient of <5(r") in the jet function J™ \ T a)- The 
sum of naive real and virtual Wilson line diagrams are integrated over the region J in the q~ , 
plane. The sum of real and virtual zero-bin subtractions are integrated over j7o, and the resulting 
sum of naive diagrams and zero-bin subtractions over the region J . Integrals over J have only UV 
divergences as long as a < 1. For a = 1, an IR divergent region remains. 



The shape of the momentum regions contributing to the one-loop soft function in Fig. 2 
suggest a simple physical interpretation of the breakdown of factorization as a — > 1. In 
the k^ plane, the region of integration in the sum of real and virtual graphs for a = 1 
is the region above the line k + k~ = Q 2 . For angularity soft functions with a < 1, as 
k^ — > oo, the loop integral goes over a region with k + k~ strictly greater than Q 2 , and in 
fact, k + k~ — > oo, while for a > 1, the loop integral enters the region with k + k~ < Q 2 , 
and in fact, k + k~ — > 0. But this latter region, k + — > oo while k~ — > or vice versa, is the 
region where collinear modes live, illustrated in Fig. 5. This means that collinear modes 
still contribute to the soft function even after the attempted factorization. 

This suggests that for a > 1, the contributions of SCETi soft and collinear modes to 
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Figure 5: Scaling of SCET modes appropriate for angularities r a , a = 0, 1. For a = 0, the collinear 
modes dominating the r a distribution have virtualities p 2 ~ (QA) 2 , parametrically separated from 
the soft scale p 2 ~ (QA 2 ) 2 . These scalings correspond to the effective theory known as SCETi . For 
a = 1, the collinear modes in the distribution have typical p 2 ~ (QA 2 ) 2 , coinciding with the soft 
scale. The collinear and soft modes are no longer separated by virtuality but instead by rapidity. 
These scalings correspond to SCETn • Collinear modes dominating angularity distributions for 
other values of a between and 1 live at scales intermediate between these limits. 

the angularity distribution have not actually been separated. In SCETi , soft, collinear, 
and hard modes can be distingushed by their well-separated virtualities, namely, p 2 s ~ 
(QA 2 ) 2 , p 2 j ~ (QA 1 ^ 1 ^/ 2 )) 2 , and p 2 H ~ Q 2 . At a = 1, the virtualities of soft and collinear 
modes contributing to the r a distribution coincide, and SCETi must be matched onto 
SCETn where collinear and soft modes both have virtualities p 2 ~ (QA 2 ) 2 . In this case, 
the modes are no longer distinguished by their virtuality, but instead by their rapidity, as 
illustrated in Fig. 5. Ref. [28] suggested a modified version of the factorization theorem 
Eq. (1.1) in which soft and jet functions are defined either with cutoffs on rapidity or in 
dimensional regularization with the scale fi separated into two light-cone scales n^, which 
must satisfy = /i 2 , with each of the two jet functions depending on one of these 

scales, and the soft function on both. However, in the present paper we do not pursue such 
a strategy and limit our analysis to angularities with strictly a < 1 . For arbitrary values of 
a, the virtuality of collinear modes p 2 ~ (QA 1 ^ 1- "/ 2 )) 2 suggests an interpretation as the 
modes of an effective theory "SCETi +a ." 5 Since our analysis and calculations utilize the 
framework of SCETi , we may expect non-negligible corrections to our results to arise for 
values of a less than but approaching 1, and for reasonable criteria for when corrections 
are negligible, our analysis is reliable for values of a < 1/2 [23]. 

4. NLL Resummation of Logarithms and Fixed-order Matching to QCD 

The fixed-order NLO cross-section, obtained by using the fixed-order expressions for the 
5 We would like to thank M. Strassler for suggesting this terminology to CL. 
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hard, jet, and soft functions in Eqs. (3.2), (3.15), and (3.35) in the factorization formula 
Eq. (2.13), contain logarithms of fi divided by the scales Q, Qr a , and the intermediate scale 
Qr a 1/(2 " a) . This means that there is no single choice for the scale [i that will simultaneously 
set all of the logarithms in the NLO cross-section to zero. For small r a , these scales become 
widely separated and the logarithms of ratios of these scales become large, which causes 
the perturbative series to break down. In Sec. 4.1 and Sec. 4.2, we take advantage of the 
effective theory framework separating the hard, jet, and soft contributions by evolving each 
of them separately through renormalization-group (RG) evolution which resums these log- 
arithms. We then combine these RG-evolved functions into the full cross-section accurate 
to NLO at fixed order in a s and resummed to NLL accuracy in Sec. 4.3. 

Since our final result for the NLL/NLO resummed distribution is derived using an 
effective theory which is valid only in the small- r a limit, it does not get the larger- r a region 
as accurately as QCD at 0{a s ). To arrive at a result that retains NLL/NLO accuracy in 
the small- r a region while retaining the accuracy of QCD at 0{a s ) in the larger- r a region, 
we need to match our distribution onto QCD. This matching is constructed such that if 
we turn off the resummation, the distributions should agree with full QCD to 0{a s ). We 
perform this matching in Sec. 4.4. 

4.1 Hard function at NLL 

The anomalous dimension of the hard function in Eq. (3.2) can be found by requiring that 
matrix elements of the bare two-jet operator in Eq. (3.3) are independent of the scale /x, 
and is given by 

7ff(A0 = -7o(Qn/2, -Qn/2; M ) - 7 o(Qn/2, -Qn/2; y) , (4.1) 

where 

lo(S>n,Pn,V) = -Zo 1 (Pn,Pnlv)V^Zo(Pn,Pn,IJ') = f 2 In — ^— ^ + 3^) , (4.2) 

d/i 27T V 2p n ■ p n J 

so that 



7^) = -^(21n^+3) , (4.3) 



which is the first term in the expansion of the anomalous dimension to all orders in a s , 

E 
Q 



7#(a0 = T H [a s ] In ^ + 1h[ol s ] ■ (4.4) 



Solving the RG equation, 



u^H(Q;u)= lH (u)H(Q;u), (4.5) 



for H(Q; u) gives 

H(Q;f,) = H(Q;no)e K » , (4.6) 
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where ojh and Kh are denned as 
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Here r 



2^1) i and 0o, 0i are the one-loop and two-loop coefficients of the beta function, 
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The two-loop running coupling a s (fi) at any scale is given by 



1 



1 



ft In f M I A 
a s (n) a s {M z ) ' 2tt n \M Z ) 4tt0 



+ 



In 



l + ^a s (M z )ln 



(4.8) 



(4.9) 



(4.10) 



In Eq. (4.6), we have used the fact that to all orders in perturbation theory, r#[a s ] is 
proportional to r cusp [a s ] , where 
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The ratio of the one-loop and two-loop coefficients of r cusp is [67] 



r 1 
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(4.11) 



(4.12) 



T,, 1 and 0\ are needed in the expressions of uijj and -fT# for complete NLL resummation 
since we formally take a 2 lnr a ~ 0(a s ). 

4.2 Jet and soft functions at NLL 

The jet and soft functions obey the RG equation 

-+oo 



/i—F(t; fj,) = / dr -y F (r - r ; /u)F(r'; ft) , 



(4.13) 



where F = J,S. The anomalous dimensions 7^ can be found from the Z-factors (given in 
Eqs. (3.16) and (3.37)) via the relation 



7 f(t - r'; /x) = - J dr" Z F \r - r"; M ) /x^^(r" - r'; /i) 



(4.14) 
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We find that 
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Both anomalous dimensions are the first terms in the perturbative expansion of the general 
form to all orders in a s [26, 68], 
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Table 1: T F , and jV for the jet 
and soft functions. 



We summarize the coefficients T F and 7^ and the j'i?-values for the jet and soft functions 
in Table 1. 

The solution of the RG equation Eq. (4.13) with the anomalous dimension tf of the 
form given in Eq. (4.17) with particular values of j F was developed in the series of papers 
[51, 69, 70, 71]. Later, it was solved for arbitrary j F in [26] using a convolution variable 
t = Q 3 r with mass dimension j = j F . The resulting evolution equation for F is 



(4.20) 



F ( T 'i M) = J dr ' Uf(t - r; [i, / u )F(r / ; /i ) , 
where the evolution kernel U F is given to all orders in a s by the expression 

U f (t - r';/x,/i ) 
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where 75; is the Euler constant and where uj f and K F are defined as 
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The plus function in Eq. (4.21) for all oj < 1 and u> ^ is defined as e 
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with the latter plus functions defined in Eq. (3.13). 

For the NLL parameters of the evolution kernel Uf, Eq. (4.22) gives 
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where we have used the fact that Tp oc r cusp . This proportionality is well known for the 
a = jet and soft functions. In Appendix A we verify that it remains true for all a < 1. 

From Eq. (4.20) we can write explicit formulas for the resummed jet and soft functions 
at any scale \x. Details of evaluating the integral over the convolution variable t' are given 
in Appendix B. For the soft function, we plug the fixed-order NLO result Eq. (3.15) at the 
scale fio into Eq. (4.20), and obtain at the scale fi, 
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and for the jet function, plug in the fixed-order NLO result Eq. (3.35) at //o into Eq. (4.20), 
and obtain at /i, 
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where in the above two equations Kp = Kp([i,fj,o), ujp = u>p(fi,fio), H(z) is the harmonic 
number function, and ip( u \z) is the polygamma function. 



6 Note that from the definition in Eq. (4.23), for uj < the ' + ' label can be dropped and so Eq. (4.23) 
is consistent with the distribution relation Eq. (3.12). 
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4.3 Full distribution at NLL 

By running the hard, jet, and soft functions from the scales /io = Hh-> MJj an d Ms> respec- 
tively, to the common factorization scale fi and performing the convolution in Eq. (2.13) 
(see Appendix B for details), we find for the final resummed expression for the two-jet 
angularity distribution with NLL/NLO perturbative accuracy 



l_do/ T 
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Q = 2 ujj{h, fij) + ujs(t-i, Ms) 

K = K H (n, fj, H ) + 2Kj(^i, fij) + Ks(p, fJ-s) , 

with u>h-,Kh given by Eq. (4.7) and ojj t g and Kj t s given by Eq. (4.24) and 
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and /(a) was defined in Eq. (3.36). 

From these expressions, it is clear that the logarithms are minimized by choosing fj,jj, 
fj,j, and (is of order Q, QT a l ^ 2 ~ a \ and Qr a , respectively. We will describe in more detail 
precisely which values we choose for these scales when we plot the full distributions in 
Sec. 6. 

4.4 Matching to QCD 

One way to achieve matching onto QCD is to include three-jet operators in the matching 
of the QCD current onto the SCET operators in Eq. (2.8) [44, 58]. For the scope of this 
paper, however, we simply adopt the matching procedure described by [42] , as implemented 
in [11]. 

To 0{a s ) the full QCD distribution will take the form 
1 da 



a dT a V27T 



a. 



A a ( Ta ) + <D(aj) . 



(4.32) 
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In Appendix C we describe how to calculate A a {T a ) numerically. Meanwhile, the fixed-order 
two-jet angularity distribution in SCET at 0{a s ) is given by the convolution Eq. (2.13) 
of the fixed-order hard, jet, and soft functions Eqs. (3.2), (3.15), and (3.35). The result is 
independent of fi (except through a s = a s (fj,)), and is given by 

5(r a )D s a + ^[D a (r a )] + , (4.33) 



a dr a 2tt 
where 

+ 4 f dx x2 ~ 2X + 2 lufx 1 " + (1 - x) 1 -"] j (4.34) 
Jo x J 

2C F g(r a )(3 + 41nr a ) 

D a {Ta) = ~- . (4.35 

2 - a T a 

The two-jet fixed-order SCET distribution Eq. (4.33) reproduces the most singular 
parts of the full QCD distribution 7 Eq. (4.32), that is, the coefficient of the <5(r a ), l/r a 
and (l/r a )lnr a pieces. The expression for D a {T a ) in Eq. (4.35) makes explicit that the 
angularities are not infrared-safe for a = 2. 

The difference of the two fixed-order distributions Eq. (4.32) and Eq. (4.33) away from 
T a = is a purely integrable function, 

By adding this remainder function to the NLL resummed SCET distribution, we obtain a 
result which both agrees with QCD to 0{a s ) and resums large logarithmic terms in the 
entire perturbative series with NLL/NLO accuracy. The matched distributions are thus 
defined as 



1 da PT 1 da 2 PT 



+ r a (r a ) . (4.37) 

NLL/NLO 



NLL/NLO °"0 dr a 

To find r a (r a ), we numerically obtain A a (T a ) from an analysis of the full QCD distributions 
away from r a = using the procedure described in Appendix C, and then subtract out the 
expression for D a (T a ) given in Eq. (4.35). 

For the case a = (thrust), the analytic form of d<r PT /dro is known [41], with which 
our formula Eq. (C.10) for ^4o( T o) agrees. Using Eqs. (C.10) and (4.35), we obtain the 
remainder function 
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which we see is integrable down to to = 0. 

7 Technically, we mean that the difference of the two distributions integrated from to e vanishes as 

€ -> 0. 
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As a consistency check of this matching technique, we calculated the total integral 8 of 
our fixed-order result, 

/ 1 da 2 PT \ 
o"total = / dr a — + 

) , ( 4 -39) 
Jo V^o dr a J 

and compared with the total inclusive cross-section, a{e + e~ — > X) = <ro(l + at 8 /n). We 
found that our results agreed to any arbitrary precision which could be achieved by our 
numerical computation. 

5. Nonperturbative Model for the Soft Function 

In this section we adapt the model for the soft function used in jet mass and thrust distribu- 
tions as constructed in [31] to work for all angularities with a < 1. This model is designed 
to describe the small- r a region where perturbation theory breaks down, while leaving the 
perturbatively reliable large- and intermediate- r a regions unaffected. The gap parameter in 
this model is designed to turn off the soft function at energies below a minimum hadronic 
threshold. Such a parameter is known to have renormalon ambiguities [31], which must 
cancel those in the perturbative soft function (which we denote in this section as S PT ) to 
yield a renormalon-free total soft function S. To ensure perturbative stability, a scheme 
is needed to explicitly enforce this cancellation order-by-order in perturbation theory. Re- 
cently, the position-mass scheme developed in Ref. [72] was used to define a renormalon-free 
gap parameter for hemisphere jet masses in Ref. [48] . This gap parameter obeys transitive 
RG evolution and has a well-behaved perturbative expansion. We implement this scheme 
generalized to arbitrary angularity. 

5.1 Review of hemisphere and thrust soft function models 

To motivate the functional form of the model function that we will use for all angularity 
distributions, we begin with the model hemisphere soft function constructed in [73]. This 
model is a function of two variables which can be chosen to be l + and l~ , defined as the 
+ and — components of the momentum in the n and n hemispheres, respectively. It takes 
the form 

The parameter A controls how steeply the soft function falls as — > 0, and B contains 
information about the cross-correlation of the soft particles in the two hemispheres. / exp is 
normalizable for A > and B > — 1. A is an O(Aqcd) parameter that describes the range 
that hadronic effects can smear the soft function around a given Finally, M(A,B) 

is chosen such that / exp is normalized to unity, dZ + d£~/ cxp (/ + , l~) = 1. 

s The upper limit on r a in Eq. (4.39), T™ ax , is that of the maximally symmetric three-jet configuration, 
r sym (a) = l/3 1 ~ a/2 [38], but only for a > -2.6 (see Appendix C). 
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In Ref. [26], this model was used to relate the total hemisphere soft function Shemi^ -1- , I ) 
to the perturbative hemisphere soft function 5'hemi(^ + ' '~) via * ne convolution 

/+oo 
dF+dr^l i(z+ - f+, r - r;n)f e *v(£ + - a, r - a) . (5.2) 
-00 

where A is the gap parameter. This method of implementing the model function ensures a 
smooth continuation between the nonperturbative, model-dominated and the perturbative 
regions of the cross-section. 

To use this expression in our formalism, we first relate the a = soft function, So(tq, //), 
and the hemisphere soft function, Sh em i(l + , l~ , fi) . Using that To = (l + + l~)/Q, we find 

So(r ; fi) = J dZ + dnShcmi(? + , l~',fJ-)S (jo ~ ^ ^ 

= Q Jdl S hcm[ (l,QT - I; fi) . (5.3) 
This gives the model function convolution for <So(to; h) as 

S (r Q ; fi) = Q jdl J dl+dl' SlJ mi (l - l+, Qro-l- r ■ /i)/ cxp (Z+ - A, T - A) 

= j dr> 5 PT (r - 4; /,)/ cxp (t[> - ^) , (5.4) 
where (absorbing A and B dependent constants into the normalization AT) 
/™"(t) = Q 2 J dr' f'-"(Qr - Qt', Qt>) 

= B)f (£ )""' Q, i + A, (B - e-<-^ . (5.5) 

/ exp (r) inherits its normalization from f exp (l + J^ dr/ Cxp (r) = 1. 
5.2 Adaptation to all angularities 

For nonzero a, we still want to use a convolution of the form 

S a (r a ; n) = J dr' a S p J(r a - r' a ; M ) / a cxp - ^ . (5.6) 

Moreover, we would like to retain the functional form of / exp since it has had relatively good 
success in describing different event shapes with the same values of A and B [73] . However, 
we must at a minimum modify / exp so that the first moment of S a {r a ; fj,) satisfies the scaling 
relation given in Eqs. (2.19) and (2.22). In terms of the first moment of 5 PT (r a ;/x) and 
/a Xp , the first moment of 5 a (r a ; n) is 

j dr a r a S a {T a ; /i) = J dr a r a j dr' a S? T {T a - T a ; fi) /^ xp (r' a 
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S P J [1] M + 2 -^ + fr [1] , (5-7) 
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where here Sa^^i^i) and /a are the first moments of (r a ; (j) and /a Xp (r a ), respec- 
tively, and in the third line we dropped a s corrections to the 0(Aqcd/<5) power corrections 
A a /Qand/ a cxp[1] . 

Since the first moment of the perturbative soft function, S a , already obeys the 
proper scaling (cf. Eq. (3.15)) we simply rescale the gap parameter, 

Aa = — - , (5.8) 
1 — a 

and require that the parameters of /a Xp vary from those in / exp such that 

/exp[l] = f dTa Ta fT V {Ta) = J_f dTT r v (r) = J_ jexpll] (5 9) 

J I- a J I- a 

This latter condition is most easily satisfied by fixing A and B to their value at a = 
and allowing A — ► A a to vary accordingly. Note from the definition of / cxp , Eq. (5.5), 
A/ cxp (Ar/Q) is independent of A and hence A a /r p (A a r/Q) = A/ cxp (Ar/<2) when A and 
B are fixed. This implies that 

rPlll = (^)7*-.r.ff- (£r.) = (^) /drrr* (£r) = (£) /«* , 

(5.10) 

and so to satisfy Eq. (5.9) we take f% xp to be defined as in Eq. (5.5) but with A replaced 
with A a where 

A 

Aa=- . 5.11 

1 — a 

5.3 Renormalon cancellation 

We want to ensure that the 1/Q renormalon ambiguity in l S' PT (r a ;/i) is cancelled order- 
by-order in perturbation theory. To implement the position-mass renormalon cancellation 
scheme defined in Ref. [72] for jet-masses and applied to the a = gap parameter in 
Ref. [48], we first take the Fourier transform of S a (T a ; fx) with respect to Qr a , 

S a (x a ; fi) = J dr a e- iQTaXa S a (r a ; M ) 

= / dr a J dr> S PT (r Q -r' a ^) / a cxp 

= S a PT (w) r xp (x a )e" 2lA ^ 

S PT (x a ; M ) e - 2 ^(MKj J / ex P(xa)e -2iA aM x a j (5 12) 

where in the second line we used Eq. (5.6) and in fourth line we split A a into two fi 
dependent pieces, A a = A a (/z) + £ a (/•*)• Note that since A a is \i- independent, S' PT and S a 
obey the same RG equation. 

Next, we demand that for some value R, the term in the first pair of brackets in the 
last line of Eq. (5.12) satisfies 

ln[S PT (* a ; ^e- 2i5 ^} iXaHRe , E) -, = , (5.13) 
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a condition which guarantees no ambiguity in S PT at order 1/Q. This gives 5 a (fx) to all 
orders in terms of S PT (T a ; fx) as 



2 JdT a e-eW(^)sPT (Ta;/x) ' 
which to leading order is given by the expression 

Since A a = A a (/i) + (5 a (/x) is /U-independent we find that to 0(a s ), 



/i— A a (/t) = <5 a (/i) = i?e 7B 

d/i d// 



1 - a V 47r 



o ( <Xs([J>) 



47T 



(5.14) 



(5.15) 



(5.16) 



Using that T^[a s ] oc r^fa^] (cf. Refs. [48, 72]) to all orders, we find that the NLL expression 
for ixdA a /d(x is and that, for arbitrary a, Ts[a s ] oc r cusp [a s ] (cf. App. A) 



M^M//) = -Re* 



r ( a s(v) \ {-, , r cusp 



«s(M) 



4tt 



which has the solution 



Re lE 

A tt (/x) = A a (^ ) —us(ji, fj, ) 



(5.17) 



(5.18) 



where //o) is given in Eq. (4.24). Note that since <5*(/x) and A a (/i) — A a (/io) are 

proportional to 1/(1 — a), Eq. (5.8) suggests that we should choose A a (fXo) to be A(/xo)/(l — 
a), where A(/xo) is the best choice for a = 0. 

Expanding Eq. (5.6) in powers of a s to 0(a s ) gives 



5 (r ;/i) = y dT^ 



^ PT (r a - <; fx) + ^J^(r a ~ r' a , fx) 



/exp 



./ 2A»(/i) 



Q 



(5.19) 



where S^ T at NLO in the first term in brackets and at LO in the second term should be 
used since 8\ is 0(a s ). Using the fixed-order expression S PT (T a ; fx) = 5(r a ) + 0(a s ) in the 
second term and integrating this term by parts gives 



S a (T a ;fx) = J dr' a 



S P a T (ra-r' a ;fx)f^ r' a 



j 2A ( M ) 



Q 



d rexp /_ 2A a (fx) 



Q dr a 



/r p r a 



Q 



(5.20) 



Evolving S a (T a ] ix s) to the scale fx with Us{r a — T' a ; /x, fxs) as in Eq. (4.20) gives 

2A, 



Sa(Ta; /*) = / d< [5 a PT (r a - r^; „) / a cxp (t£ - 



q <Js{Ta ~ T a ] fX, Us) -foJJa [ T a 



, 2A a (in 



Q 



(5.21) 
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Here we keep A a and 5 a at the scale (is which is needed to achieve the 1/Q renormalon 
cancellation [74]. 

Finally, Eq. (5.21) implies that the total resummed distribution at NLL convoluted 
with the model function f% xp is 



1 da 

0"o dr a 



NLL/NLO 



1 da 



PT 



dT 'a\ — J- ( T a-r' a ,fl) 
I (To QTa 



fexp I I 
J a [ a 



Q 



NLL/NLO 
U a( T a ~ T' a ;Li,VH,fJ<J,fJ<s) 



2A a (/x 5 ) 



Q 



dr. 



I J a 



cxp 



2A a ((i s ] 



Q 



(5.22) 



where the resummed two-jet distribution matched to QCD, dcr^^/dT a | nll /nlo ; is given in 
Eq. (4.37) and U% is given in Eq. (4.28). 

5.4 Numerical results for the soft function 

By plugging the partonic soft function Eq. (3.15) into the model Eq. (5.20), we obtain for 
the full convoluted model soft function to 0(a s ), 



1 



+ 



a s Ci? 1 
2vr I- a 

2a s C F 1 



In 2 



7T" 

~6~ 



/r p (^ A )-^^/ cxp (^ A ) 



it l-a./o t> VQV 



[f C a 



cxp ^ A -r')-/r p (^)] , 



(5.23) 



where = r a — 2A a (//)/Q. To integrate against the plus distributions in Eq. (3.15), we 
used the prescription 



dr 



8{x) 



dx 



0{x) lnx 



9(x) 



[/(x)-/(0)] + /(0)lna 



dx- 

x 

l dx d -^^[f(x)-f(0)] + bm^a, 
o x 1 



(5.24a) 
(5.24b) 



which correspond to the definition of plus-functions given in Eq. (3.13). To minimize the 
logarithms in the peak region of the soft function while also avoiding the Landau pole in 
a s , it is natural to choose the scale to be of order /x > Aqcd- To minimize the logarithms 
for larger values of r a , it is natural to choose \i ~ Qr a . A scale choice that interpolates 
between these two regions is 



vf n )(QT a -vg m ) 2 + (vf n y 



(5.25) 



where the minimum scale is /i™ 111 > Aqcd- 

In Fig. 6, we plot S a (r a ; /j.) for six values of a between —2 and 1/2. In each plot, we show 
the tree-level (LO) soft function with a gap parameter (solid gray), the one-loop (NLO) 
soft function with a gap parameter but without renormalon subtraction (dashed green), 
and the one- loop soft function with a gap and renormalon subtraction (solid blue) . For the 
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Figure 6: Angularity soft functions with a gap parameter, at tree-level (solid gray) and at one-loop 
with (solid blue) and without (dashed green) renormalon subtraction, for Q = 100 GeV, for several 
values of a as labeled on each plot. The variation of the soft functions with the scale [i is illustrated 
by first setting /i™ 111 = 1.0 GeV in Eq. (5.25) and choosing /i to be (0.8, 1, 1.2) times the formula in 
Eq. (5.25), with the plots for smaller values of \i$ peaking earlier in r a . For the model parameters 
we take A = 2.5, B = —0.4, A = 0.55 GeV. In the renormalon subtraction Eq. (5.15), we have 
chosen R = 200 MeV. 

parameters in the model function Eq. (5.5) we take A = 2.5, B = —0.4, A = 0.55 GeV, as 
extracted from a fit to the jet mass distribution [73]. For the scale dependence of the gap 
parameter, we choose Ao(l GeV) = 100 MeV and use Eq. (5.18) to evolve to other scales. 
We choose R = 200 MeV in the renormalon subtraction Eq. (5.15) and the minimum value 
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of the scale in Eq. (5.25) to be ^t™ 111 = 1 GeV. We illustrate the variation of S a (T a ; fi) with 
the scale \x by varying it between 0.8 and 1.2 times the formula in Eq. (5.25). The tree-level 
soft functions depend on [i only through the gap parameter A a (/i) and thus artificially have 
smaller scale variation than the one-loop soft functions, at which order the nontrivial \x 
dependence is first probed. 

The one-loop soft functions in Fig. 6 display unphysical behavior near r a = by taking 
negative values, due to the renormalon ambiguity in the perturbative series for the partonic 
soft function. By cancelling the renormalon ambiguity between the partonic soft function 
and the nonperturbative gap parameter A a through Eq. (5.20), we obtain the renormalon- 
free one-loop soft functions. One of the plots of the soft function for a = 1/2 still exhibits 
a small negative dip after renormalon subtraction, but it is nevertheless much smaller than 
the original negative dip, and from its size may be expected to an effect of higher-order 
power corrections. The dip does not appear in the total cross-section calculated below in 
Sec. 6. 

6. Numerical Results for the Full Distribution 

In this section we plot the angularity distributions da/dr a which include LO and NLO 
perturbative hard, jet, and soft function contributions, resummation of large logarithmic 
terms to NLL accuracy, matching to QCD at 0(a s ), and the effects of the nonperturbative 
gapped soft functions. 

In Fig. 7 we plot the angularity distributions given by Eq. (5.22), plugging in the NLL 
resummed partonic distribution given by Eq. (4.27) and matched according to Eq. (4.37). 
We keep the same soft model function parameters as in the previous section. As noted 
earlier, the logarithms in the hard, jet, and soft functions are minimized by choosing 

= Q, fJ-j ~ QT~a^ 2 , and us ~ Qi~ a . In order to avoid the Landau pole in a s as 
T a — > 0, we choose the scales as in Eq. (5.25) , 

fi s = ^0(Qr a - nf*)(Qr a - ^ in ) 2 + (^ in ) 2 (6.1a) 
ix j = ^(Qr a 1/(2 " a) - /xf n )(Qr a 1/(2 - a) - + (/^ . (fUb) 

We may vary independently, or choose them in a correlated fashion suggested by their 
natural scaling fxg ~ Q^,Hj ~ Q\ l ^ 2 ~ a \ that is, 

^™ in = Q(l-«)/(2-*)( M ™*)V(2-a) _ ( 6 2 ) 

In Fig. 7 we have done the latter. The NLL /NLO distributions exhibit negative values for 
small Tq, clS cl result of the renormalon ambiguity. Performing the renormalon subtraction 
in the soft function removes this pathology. 
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Figure 7: Angularity distributions at Q — 100 GeV for six values of a between —2 and 1/2. The 
solid gray curves are the LO partonic distributions resummed to NLL and convoluted with the 
gapped soft model function. The dotted green curves are NLL/NLO convoluted with the gapped 
soft function but without renormalon subtraction. The dashed red curves are the same as the 
green but with renormalon subtraction, and the solid blue curves are the same as the red but 
matched to fixed-order QCD at 0(a s ). We choose the scales [i = Q,^™" 1 — 1 GeV, and /i™ m 
given by Eq. (6.2). For the gap parameter we take Ao(l GeV) = 100 MeV and in the renormalon 
subtraction R = 200 MeV. 
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Figure 8: Angularity distributions at Q — 100 GeV. The full, NLL/NLO resummed, renormalon- 
subtracted distributions in Fig. 7 are here shown all on the same scale. The parameters are chosen 
the same as in Fig. 7. From highest to lowest peak value, the curves are for a = —2, —1, — 1,0, ^, i. 
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Figure 9: Hard scale variation (dark green band) and correlated jet and soft scale variation 
(light blue band) of the NLL/NLO resummed, renormalon-subtracted angularity distributions at 
Q = 100 GeV for a = —1, a = 0, a = 1/4, and a — 1/2. For the hard scale variation, [i^ varied 
between Q/2 and 2Q and for the correlated scale variation, \ij and [is ar e varied between half the 
values given in Eq. (6.1) and twice these values. 
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In Fig. 8 we plot angularity distributions for the values of a used in Fig. 7 on the 
same figure to illustrate clearly how they change with a. The range of r a populated by 
two-jet-like events grows with increasing a, so that the peak regions are populated by jets 
of increasing narrowness with increasing a. This is reflected in the scales fij t s in Eq. (6.1) 
drawing closer as a grows to 1. 

In Fig. 9 we vary the hard, jet, and soft scales and plot the resulting variation of our 
final predictions for the distributions. First we vary the hard scale \xn between Q/2 and 
2Q, plotting the result in the dark green band. Then we vary the collinear and soft scales 
Hj : s between half and twice the values we chose in Eq. (6.1) and plot the result in the light 
blue band. 

Although published data on e + e~ angularity distributions for a ^ are not yet avail- 
able, data for the a = (thrust) distribution are of course plentiful. The remaining differ- 
ence between our prediction in Fig. 7 and existing measurements of the a = distribution 
can be accounted for by higher-order perturbative corrections (see, for example, Fig. 6 in 
Ref. [11]), which are known but have not been included here, since we calculated the other 
angularity distributions only to NLL/NLO. For a sufficiently smaller than 1, we expect our 
predictions of all angularity distributions to agree with data to the same accuracy that the 
NLL/NLO a = prediction agrees with the thrust data. 

7. Comparison to Previous Results and Classic Resummation 

To compare to previous predictions of angularity distributions [10, 38] and focus more 
generally on the differences between SCET and alternative approaches to factorization and 
resummation, in this section we restrict our attention to the perturbative distribution both 
before matching, Eq. (4.27), and after matching, Eq. (4.37), leaving out the nonperturbative 
model of Sec. 5. 

Our result for the unmatched NLL resummed distribution Eq. (4.27) involves an evo- 
lution factor Uq, which resums all leading and next-to- leading logarithms (for example the 
(1/T a )lnr a and l/r a terms in the fixed-order D a (r) of Eq. (4.35)), and a multiplicative 
NLO prefactor 1 + f H + 2/j + f s = 1 + 0{a s ). Both the evolution factor and the NLO 
prefactor are sensitive to physics at the three distinct scales fJ-j, and fis- Keeping these 
scales arbitrary until after solving the RG equations in Sec. 4 and retaining the freedom 
to choose them only at the end provides a flexibility which is indispensable in achieving 
reliable predictions in the SCET approach. This approach has significant advantages over 
what we refer to as the classic approach to resummation in QCD [42] . 

To illustrate these advantages, we compare our results for angularity distributions to 
those obtained in full QCD [10, 38]. The analysis in Ref. [10] used a formalism of factoriza- 
tion and resummation of logarithms through renormalization-group evolution paralleling 
that of SCET, in principle containing all the advantages that we emphasize here, but which 
were not fully realized. Before arriving at the explicit prediction for the NLL resummed 
distribution dcr/dr a given in Ref. [38], the factorized result of Ref. [10] was first converted 
into the form of a resummed event shape distribution that would be obtained using the 
classic approach (and has been for a = 0). 
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Figure 10: Factorization scale /i variation of the (unmatched, partonic) SCET NLL/LO (light 
blue band) and the classic QCD NLL/LO (red band) resummed results for angularity distributions, 
/i is varied over the range ^ < /i < 2Q with Q = 100 GeV for the cases a = —1, a = 0, a = 1/4, 
and a = 1/2. To make a direct comparison to the QCD results, the scales in the SCET results have 
been chosen as /i = /i# = Q, H,i = QtV^ 2 , and /is = Qr a . 



One major advantage of the SCET approach over the classic approach is the presence 
of Landau pole singularities in the results of the classic approach that are not in the re- 
sults from SCET, as also found in the cases of DIS and Drell-Yan [49, 50, 51]. We can 
illustrate why SCET avoids this for the case of angularities by returning to our results for 
the resummed jet and soft functions and for the final resummed distribution. From the 
expressions for the resummed soft function S a (r^), Eq. (4.25), and for the resummed jet 
function J™(r"), Eq. (4.26), one might be tempted to set y, s = Qt^ and \ij = Q(To) 1 ^ 2-0 ^ 
since the logarithms in Eqs. (4.25) and (4.26) are minimized for these choices. The prob- 
lem with this choice is that the soft and jet functions still enter the convolution in the 
factorization theorem Eq. (1.1) and thus the scales in a s {^j/s) run below Tq' s = Aqcd/Q 
even for r a > Aqqd / Q (where r a = r™ + + ) if these r" ' s -dependent scales are chosen. 
However, for a r^-independent choice of /is in the case of the soft function, for instance, 
the full functional dependence of the resummed S(t^; fJ,) on r* and /is is such that after the 
integrals over r*, r™, and needed to get to the final resummed distribution, Eq. (4.27), 
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are performed, the resulting dependence on us only comes in the combination ^s/Q T a m 
logarithms (and similarly for the jet functions). The proper choice is thus fis ~ Qi~ a (and 
/xj ~ Qt\^ 2 and not [is ~ Qt^. With this choice, Landau pole singularites never affect 
our result for r a > Aqcd/Q- Setting [is = Qt^ before doing the convolution Eq. (1.1) is 
equivalent to setting [is = Q/v in the Laplace transform with respect to v of the distri- 
bution, which is the scale choice made in Ref. [10] needed to reproduce the classic result 
for a = 0. Thus, when transforming back to get dcj/dr a , one inevitably runs into spurious 
Landau pole singularities with this scale choice 9 , confirming the similar observation of [51]. 

Another difference between the explicit results we give and those given in [10] is that 
while both achieved resummation of logarithms to NLL accuracy, the latter does not include 
a full NLO calculation of the jet and soft functions in the distribution da/dr a , that is, 
effectively does not have the prefactors Jh.j,s- As with our SCET results, the results of [10] 
are not as accurate as fixed-order QCD in the large- r a region and need to be matched. This 
matching was subsequently performed numerically to at 0(a 2 s ) in Ref. [38]. We summarize 
this by saying that we have resummed logarithms of r a to NLL/NLO with 0{a s ) matching 
and Ref. [38] has resummed to NLL/LO with 0(a 2 s ) matching. 

The explicit dependence of the NLO prefactor on the separate scales [Ih,j,s makes it 
distinct from what is obtained by NLO matching to QCD in the large- r a region where 
the three scales are comparable. Specifically, it improves the accuracy in the smaller- 
T a region where the distribution depends on physics at the three widely disparate scales 
separately, as revealed by the factorization theorem. We emphasize that even though 
the effects of including this NLO piece are formally of next-to-next-to-leading logarithmic 
(NNLL) accuracy (using the counting a s lnr a ~ 0(1)), it is natural to include it in our 
NLL resummed result since the dependence on the arbitrary scales [Ih,j,s is cancelled to 
order a s in our NLL/NLO calculation. 10 

Finally, we point out that while SCET can incorporate 0(a 2 ) matching with, for 
example, an O(o 2 ) QCD calculation or an event generator, the classic approach by itself is 
less easily generalized to achieve full NLL/NLO accuracy. The reason for this difference is 
that SCET predicts the evolution boundary conditions for the hard, jet, and soft functions, 
H(Q; /if/) and -F(r a ; fip) (F = J, S) in Eq. (4.20), for arbitrary scales [Ih,j,s order by order 
in perturbation theory. On the other hand, as discussed in Ref. [45], the classic approach 
in contrast must effectively use the evolution boundary conditions F(r a ; /io) = 5(T a ), which 
are LO in the SCET point of view. An implication of this difference is that, since our NLO 
prefactor is formally part of the NNLL series, full NNLL resummation is a nontrivial task 
in the classic approach (e.g. [76, 77]) whereas it is straightforward in SCET, using no new 
techniques additional to the ones described above. 

9 There are also inherent Landau pole singularities in the classic approach before transforming back to 
T a -space and thus not associated with making ^-dependent scale choices for fij t s- ln the classic approach, a 
prescription to avoid both types of Landau pole singularities is employed, but at the expense of introducing 
unphysical power corrections [42, 75]. The results of [38] plotted in Fig. 10 used the prescription of [42]. 

10 More generally, in an N Tl LL/N m LO calculation, the dependence on fiH,j,s cancels up to order Q ™ ln { n ' m } ; 
as the [ih,j,s derivative of the logarithm of the distribution receives contributions from the prefactor at 
order a™ and from the anomalous dimension at order a". 
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In Fig. 10, we compare our result with the classic result obtained in [10]. To make this 
comparison, we truncate our result to NLL/LO accuracy and make the scale choices that 
are equivalent to those that were made in Ref. [10] for the purpose of arriving at the classic 
resummed form. Namely, we run the jet and soft functions from their respective natural 
scales, fij = Qt\^ 2 ^ and fis = QT a , to the hard scale set to fin = Q- In addition, in 
Ref. [10] the factorization scale \x was also chosen to be /i = effectively turning off 
running between \xn and [i. Thus, to make a genuine comparison, we vary /i both in the 
classic result given in [38] and in our result Eq. (4.27) over the range Q/2 to 2Q, fixing 
/i h = /x in our result. Notice from the plots that the peak position appears to be more 
stable in the SCET results relative to the classic results and that there is a discrepancy 
in the overall normalization in the peak region, both of which may be attributed to power 
corrections arising from the spurious Landau poles present in the classic result. 

8. Conclusions 

We have calculated angularity distributions in e + e~ collisions for a < 1 to 0{a s ) in fixed- 
order accuracy, resummed leading and next-to-leading large logarithms in the perturbative 
series, incorporated the effects of a nonperturbative model for the soft function with a gap 
parameter, and cancelled the leading renormalon ambiguities in the perturbative expansion 
of the distribution and the gap parameter. Our new results for the one-loop jet and 
soft functions for all a < 1 and the NLL resummation of logarithms of r a with explicit 
analytical dependence on the scales Hh,j,s made possible what we believe are the most 
precise predictions of angularity distributions to date. 

These predictions, especially after extension to higher orders in perturbation theory 
and resummation of logarithms, can prove useful in improving extraction of the strong 
coupling a s or the parameters of nonperturbative models for the soft function. At the 
present time, in the absence of a new linear collider, such extractions would require the 
re-analysis of LEP data to extract the angularity distributions. 

We also gain insight into the steps that will be required to predict jet observables in 
hadronic collisions, a broad range of which have been studied in [78, 79] using the classic 
approach. An SCET-based framework to factorize jet observables in this environment was 
developed in [59] . Our analysis of angularities suggests that the study of any set of jet 
observables which vary in their sensitivity to narrower or wider jets or which depend on a 
jet algorithm picking out narrower or wider jets should be scrutinized in the same way as 
we did for angularities to determine whether the contributions of collinear and soft modes 
to each observable can be clearly separated. Also, our calculations of light quark angularity 
distributions in e + e~ collisions can be extended to calculating individual jet shapes for jets 
of various origins to higher accuracy, contributing to strategies to use such jet shapes to 
distinguish experimentally different types of jets [2, 3]. 

While we have used SCET to calculate and explore the behavior of angularity distri- 
butions, the variation in behavior of the angularities has in turn shed light on the behavior 
and applicability of the effective theory. Varying a essentially varies the collinear scale of 
SCET, in effect interpolating between (and extrapolating beyond) SCETj and SCETn , 
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and so angularities provide an ideal testing ground for the behavior of these effective the- 
ories. 

It is natural and straightforward to consider further improvement of our predictions 
to higher perturbative accuracy and reduced nonperturbative uncertainty. We believe by 
using the cut diagram methods described above to obtain the angularity distributions to 
0(a s ) we can extend our results to O(a^) in a straightforward manner. Also, all of the 
ingredients necessary for NNLL resummation at a = are already known [11], and we 
would only need to calculate those pieces which change with a. The three-loop Tj t s part 
of the jet and soft anomalous dimensions for arbitrary a can be obtained from the known 
three-loop r cusp [51] and the all-orders proportionality Tj t s oc r cusp which we verified in 
Appendix A. The only unknown ingredients are the two-loop non-cusp part of the jet and 
soft anomalous dimensions. These can be obtained solely from the UV divergences of the 
two- loop graphs, and would immediately extend our results to NNLL accuracy. As for 
nonperturbative effects in the angularity distributions, we have treated these effects in the 
soft function in the simplest manner possible, adapting the a = soft model function to 
all a by rescaling its first moment. Comparison of these predictions to e + e~ data can shed 
light on the reliability of this choice. 

Angularities and other event shapes have proven to be powerful probes of QCD and its 
effective theories, and promise to play a key role in the new era of collider physics searching 
for signals of new physics amid a sea of jets and strong interactions. 

Acknowledgments 

We are particularly grateful to C. Bauer, G. Sterman, and I. Stewart for a detailed review 
of an early draft and much valuable feedback and insight. We also thank M. Fickinger, 
S. Fleming, and A. Jain for enlightening discussions. CL and GO are grateful to the 
Institute for Nuclear Theory at the University of Washington for its hospitality during a 
portion of this work. AH is supported in part by an LHC Theory Initiative Graduate 
Fellowship. This work was supported in part by the U.S. Department of Energy under 
Contract DE-AC02-05CH11231, and in part by the National Science Foundation under 
grant PHY-0457315. 

A. Relation Among Hard, Jet, Soft, and Cusp Anomalous Dimensions 

In Eq. (4.24) we used that the r^[a s ] part of the jet or soft function anomalous dimension, 
defined in Eq. (4.17), is proportional to the cusp anomalous dimension r cusp to all orders 
in a s . This fact is well known for the standard a = jet function and soft functions. In 
this section we verify that this relation remains true for all a. Our strategy will be to show 
that r j,s[a s ] must always remain proportional to r# [«<;], which is independent of a and is 
already known to be proportional to r cusp . 

The consistency of the factorization theorem Eq. (2.13) requires a relation among the 
hard, jet, and soft function renormalization counterterms, and, thus, among the anomalous 
dimensions (see, e.g., [10, 26]). This relation can be derived by requiring that Eq. (2.13) 
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remain true when written in terms of either the bare or renormalized hard, jet, and soft 
functions on the right-hand side. This requires that 

z h 1 (/ j >) s (tj ~ T s) = J dr' J &t"Zj(tj - t'; \i)Zj(t' - r"; l±)Z s {r" - t s ; fJ,) , (A.l) 

to all orders in a s . To 0(a s ), we can easily verify this relation using Eqs. (3.4), (3.16), and 
(3.37) with Zjj(fi) = \Zq{^l)\~ 2 . This relation amongst the counterterms requires in turn 
that the anomalous dimensions satisfy 

- 1h W{t) = 2 7 j(r; /x) + 7s (r; //) . (A.2) 

To all orders in a s the hard anomalous dimension takes the form of Eq. (4.4) and the 
jet and soft anomalous dimensions take the general form of Eq. (4.17) [68] , where the 
constant jp is jj = 1/(2 — a) for the jet function and jg = 1 for the soft function. The 
constraint Eq. (A.2) then requires the three independent relations 



= -Tj[a s ] + — T s [a s ] , 
3J Js 

-T H [a s ] = 2Tj[a s ] +T s [a s ] , 
-1h[ol s ] = 2-fj[a s ) + 75 [a s ] , 



(A.3) 

(A.4) 
(A.5) 



to all orders in a s . These relations can be verified to 0{a s ) from Eq. (4.3) and Table 1. 
The first two relations Eqs. (A.3) and (A.4) taken together imply that 



1 



-T H [a s ], Tj[a s ] 



1 - a/2 



(A.6) 



1 — a 1 — a 

to all orders in a s and for all a < 1. Since r#[a,s] oc r cusp and is independent of a, both 
T s ,j[a s ] oc r cusp as well. 

B. Evaluation of Resummed Jet and Soft Functions and Full Distribution 

To evaluate the resummed jet and soft functions, we used the following method. First, note 
that from the expressions for the evolution equation, Eq. (4.20), the form of the evolution 
kernel, Eq. (4.21), and the generic form of the NLO jet and soft functions, 



F(t;hq) = cxS(t) + c 2 



+ C 3 



lnr 
r 



the resummed jet and soft functions are proportional to 



/dr' 


" 6{t-t') ' 


_(t - r') 1+w _ 





F(r'; fio) = aWi + c 2 W 2 + c 3 W 3 , 



(B.l) 



(B.2) 



where 



Wi 

w 3 



= fa, 


' 9{t-t') ' 
_(t-t') 1 +"_ 


+ 






' 6{t-t') ' 
_(t - T ') 1 +^_ 
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\6(T>y 


+ 




' 6(t-t') ' 




"#(T')ln(T') 




+ 


T 





(B.3) 
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Next, note that from the definitions of the plus functions, Eqs. (3.13) and (4.23), we can 
find Wi as the coefficient of 5 % in the Taylor series of W(S), where W{5) is defined as 



W(S) = / dr 



' 9{t-t') 
{t-t'Y+^ 

Eq. (B.4) follows from the fact that 



0(r') 

t'1+5 



9{t - t" 



\T — r 



/Al+u>i 



T 



T 



r(-u)F(-6) 
T(-uj - 5) 

r(-wi)r(-u;2) 



-I+uj+5 



-U>1 — L02) 



9{t - T 1 



(T — T 



(B.4) 



(B.5) 



By expanding both sides of Eq. (B.4) in 5 and comparing like powers of 5, we find that 



Wi 

w 3 



0(r) 



w 2 



ln(r) 



ln 2 (r) -ln(r) J ff(-l - u) + 



H(-l 



7T" 

12 



-1+0) 



0(r) 



(B.6) 



Here, H(z) is the harmonic number function and ip( u \z) is the polygamma function. 

The same technique can be used to analytically calculate the fully resummed cross- 
section, Eq. (2.13), directly from the unresummed jet and soft functions. The resummed 
cross-section is of the form 



1 da 
o-q dr 



PT 3 ( f 

«]] / dTid7(Fi(T-;fM 



( r ._^)l-K* 



8{t -TX-T2- T 3 ) . 



(B.7) 



where the jet and soft functions ^(r^/ij) are all of the form given in Eq. (B.l). These 
integrals can be done most easily by replacing the Fifa; fii) on the right-hand side of 
Eq. (B.7) with [0(r)/r 1+<5i ] , expanding in <5j before and after combining all the plus 
distributions using Eq. (B.5), and comparing like powers of the 5i. The result for the 
resummed cross-section Eq. (4.27) then follows. 



C. Angularity Distribution in QCD to 0(a s ) 

In Sec. 6 we matched the NLL resummed two-jet angularity distributions in SCET onto 
the 0{a s ) fixed-order distributions in full QCD using the remainder function r a (r a ), de- 
fined in Eq. (4.36). In this section we provide some details of how we calculate the QCD 
contribution to r a (r a ) away from r a = 0, A a (T a ). In the process, we show that for a < — 1.9 
the angularities of events with more two-jet like kinematics become degenerate with those 
of more three-jet like events and contribute to the same r , and that for a < —2.6 the 
maximally symmetric three-jet event contributes to a smaller r a then some more two-jet 
like events. Thus, for small enough a, angularities fail to separate two-jet and three-jet like 
events. 

Both the one loop qq and tree-level qqg final states contribute to da/dr a at 0(a s ). 
However, the qq final states' contribution is proportional to 5(r a ) and hence only contributes 



- 44 - 



0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

Xl X\ 

(A) (B) 

Figure 11: (A) Phase space for three-particle qqg final state. The energy fractions Xi = 2Ei/Q 
of the three particles satisfy x\ + X2 + x% = 2. In region I, x\ > £2,3, in region II, xi > 2:1.3, and 
in region III, x% > xi 2- The thrust axis is in the direction of the particle with the largest energy. 
(B) Contours of constant r a = 1/10 for a = — 1 (purple), a — (gray), and a = 1 (pink). The 
differential cross-section da/dr a is given by integrals over these contours in the X1.2 phase space. 

to A s a . Thus to find A a (r a ) we only need to consider the tree-level qqg final states. Their 
contribution can be writtten as 

1 (in-™ fn \ 

where 

dxidx 27 - 1 rrz 2 r5(r a - T a (xx,X 2 )) , (C.2) 

(1 - Xi)(l - x 2 ) 

and where jci^ = 2E\^/Q are the energy fractions of any two of the three final-state 
partons. By momentum conservation, x\ + x 2 + £3 = 2. For a three-particle final state, 
the thrust axis is given by the direction of the particle with the largest energy. The x\^ 
phase space can be divided into three regions, as illustrated in Fig. 11 A, according to which 
parton has the largest energy. In the region in which is larger than Xj the angularity 
T a (xi,x 2 ) is given by 

T (X1, X 2 ) = -(1 - X,) l ~ a / 2 \{l - X] ) l - a l\l - X k T' 2 + (1 - X 3 ) a ' 2 {l - X k ) l - a l 2 . 

(C.3) 

At each fixed value of r a = c in the distribution Eq. (C.l), the delta function restricts 
the integral over x\, 2 to a linear contour determined by the equation T a {x\,X2) = c, where 
T a (xi,x 2 ) is given by Eq. (C.3). Examples of these integration contours are shown in 
Fig. 11B. 
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Figure 12: The local minimum (green line) and maximum (red line) of the function F a (w) over 
the range < w < 1/2 coincide at the point a = ai w —1.978. At a e a 2 ~ —2.618, the value of 
angularity for the maximally symmetric three-jet case, r sym (a) = l/3 1_a / 2 (blue line), intersects 
the local maximum and so for a < a2, the value of maximum angularity for such a corresponds not 
to the maximally symmetric case but to a more two-jet like event. 

It is sufficient to consider the part of the phase space corresponding to region III shown 
in Fig. If, where 23 > x\ )2 . Integration over the remaining two regions can be related to 
the integration over region III by a trivial shift of variables of integration. Thus we need 
to solve 

C = 1 (X! + X 2 - l) 1 ^'/ 2 [(1 - Xl ) l ~ a l 2 {l - X 2 ) a/2 + (1 " Xl)°/ 2 (1 - X 2 ) 1 ^ 

2 — x\ — x 2 L 

(C4) 

where x\^ 2 lie in region III. To find an explicit one-variable parameterization for x\^ 2 {w) 
which satisfies Eq. (C.4), we first absorb the factor 1/(2 — xi — x 2 ) inside the brackets and 
define 

1 — X\ 



w 



2 — x\ — x 2 
In terms of w, Eq. (C.4) can be written as 

,^/2 {1 _ w) a/2 +w a/2 {1 _ w) l-a/2 



c = { Xl + x 2 - l) 1 ^ 2 
Solving Eqs. (C.5, C.6) for xi,X2 gives: 

xi(w) = 1 — w + w 



(C.5) 



(C.6) 



w l-a/2(l _ w )a/2 + w a/2^ _ w )l-a/2 



l-a/2 



x 2 (w) = Xl(l — w). 



(C.7) 



Clearly from Eq. (C.5), w lies in the interval < w < 1. The precise range of 
values for w is determined from the conditions x\(w) < 2 — xi(w) — x 2 (w) and x 2 (w) < 
2 — x\(w) — x 2 (w). These inequalities can be simplified to 



c < mm{F a (w),F a (l - w)} 



F a (w) for 0<w<l/2 
F a (l-w) for 1/2<tij<1 



(C.8) 
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Figure 13: (A), (B), (C) Allowed regions for the parameter w as a function of fixed r a — c are 
bounded by the curves F a (w) and F a (l — w). For (A), (D) a = — 1, the integration is over a single, 
continuous domain for all fixed r a — c but for (B), (E) a = —2.3 and (C), (F) a = —4, there are 
multiple disjoint regions of integration for large enough values of c. In (D), (E), and (F), the blue, 
red, and green curves represent contours of integration for fixed r a = c, in order of increasing c, 
and correspond to integration over a range of w given by the lines of constant r a = c in the regions 
of the same color in (A), (B) and (C), respectively. 



where 

_ 7/1 W 2 

f «w - + Q- - w ^ ■ ^ 

The function F a (w) is monotonically increasing over the range < w < 1/2 only for 
2 > a > a\ ~ —1.978, but for a < a\ turns out to have exactly one local maximum, T max (a), 
and one local minimum, r m ; n (a). At a = 02 ~ —2.618, r max (a) is equal to the angularity of 
the symmetric three-jet configuration x% = %2 = X3 (where w = 1/2), T sym (a) = l/3 1_a / 2 . 
Thus, the global maximum of r a over the whole range < w < 1, defined as T™ ax , is 
7"max(o) for a < d2 and is r sym (a) for a > 02- 

In Fig. 12, we show how the maximum and minimum of the function F a (w) depend 
on a, along with the a dependence of the symmetric three-jet configuration, and plot the 
special points a\ and 02- 

In Fig. 13 we plot the boundary of r a (F a (w) for < w < 1/2 and F a {\ — w) for 
1/2 < w < 1) together with the contours of constant T a (x\,X2) = c for different values of 
c in the full xi-x? plane for the cases a = —1, a = —2.3, and a = —4, which qualitatively 
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represent the three cases a > a%, a\ > a > 0,2, and 02 > a, respectively. From this analysis 
we conclude that for a < a\ and especially a < 0,2 angularities fail to separate two-jet like 
and three-jet like events. 

To obtain A a (r a ), we evaluate the integral in Eq. (C.l) over the appropriate contours in 
the x\ 2 phase space numerically, except for a = 0, for which the integral can be evaluated 
analytically, giving (cf. [41]) 



Aq(t ) = C F 



2(2 - 3r + 3r 2 ) ^ (1 - 2t q \ 3(1 - 3r ) (1 + r 



To(l - To) V T / T 



(CIO) 
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